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SUMMARY. 

The principal difficulty in developing a dislocation theory of mechanical 
twinning has been to explain how a twin grows through a finite thickness 
of crystal. By extending the recent theory of slip bands, proposed by 
Frank and Read, it is possible to show that a dislocation can move steadily 
from plane to plane in a crystal. A theory of mechanical twinning, which 
is formally analogous to Frank’s theory of crystal growth, can then be 
developed. The theory is applied to the body-centred and face-centred 
cubic lattices. 


§ 1. INTRODUCTION. 

THE idea that mechanical twinning takes place in crystals by the contin- 
uous growth on an atomic scale of twinned material, and not by the 
simultaneous movement of all the atoms involved, has often been 
suggested ; for example, by Frenkel and Kontorova (1939), Seitz and 
Read (1941), and Frank and van der Merwe (1949). The general argu- 
ments for a dislocation mechanism of mechanical twinning are the same 
as those for the dislocation theory of slip. First, it is scarcely believable 
that the atoms concerned should all move simultaneously, and second, 
twinning takes place at stresses far below the theoretical shear strength of 
a perfect lattice. 

In some examples of mechanical twinning the atomic movements 
required to produce the twinned lattice from the original one are more 
complicated than those that produce slip. The movements can usually 
be described in terms of a homogeneous shear of all the atoms involved, 
together with localized rearrangements in which neighbouring atoms move 
in different directions from one another. The twinning shear of one plane 
of atoms over another can be produced by the movement in the twinning 
plane of a dislocation of the usual type, but with a Burgers vector that is not 
a lattice vector (an ‘“‘ imperfect dislocation ”’). On the other hand, the 
localized rearrangements cannot be described formally in terms of a 
dislocation. In practice, of course, these rearrangements will usually be 
produced automatically by the passage of the dislocation, because the 
arrangement of atoms is almost certainly mechanically unstable until the 
transformation is completed. In this paper a dislocation mechanism is 
suggested for the shear component of the twinning movements, and the 
extension of the theory to those cases where localized rearrangements also 
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take place is not considered. The theory is thus most directly applicable 
to those twinning deformations, such as that of «-iron, or shear trans- 
formations where the atom movements are completely resolvable into 
homogeneous shears on crystallographic planes. 

It is usually not difficult to discover a suitable imperfect dislocation 
to cause the required shear of neighbouring planes as it passes between 
them. The important problem is to explain how twinning develops 
homogeneously through successive lattice planes, in contrast to the 
behaviour in slip, where the deformation is concentrated on a single plane. 
In terms of dislocations the homogeneous twinning shear requires either an 
avalanche of dislocations, one on every lattice plane without exception, or 
the motion of a single dislocation successively from one plane to the next 
in aregular manner. It is the purpose of this paper to show that the latter 
process can occur naturally in certain crystals containing dislocations. 


§2. THe Pott MECHANISM. 


Frank and Read (1950) have recently proposed a mechanism whereby 
unlimited amounts of slip can be produced by a single dislocation moving 
in asingle slip plane. The essential idea is that the line of this dislocation 
does not lie for the whole of its length in the slip plane ; at some point, 
or at two points if the line does not extend to the surface of the crystal, the 
line turns out of the plane along some intersecting direction. Thus 
consider a dislocation line of arc element ds and Burgers vector b which, 
for part of its length, lies in a slip plane of unit normal v, where b. v=—0. 
Then repeated rotation in the slip plane of this part of the line about the 
interior point (or points) of emergence of the line from the plane can cause 
an amount of slip on the plane which is limited only by boundary condi- 


tions. This process requires the part of the dislocation line that leads out 


of the slip plane (which we shall call the ** pole ”’ dislocation) to be in some 
way anchored ; the moving part of the line (“‘ sweeping ”’ dislocation) then 
forms a plane spiral about the fixed point of emergence. 

It is not essential that the pole dislocation should end in the slip plane. 
Consider for example three straight dislocation lines lying along s,, s5, s3 
and meeting at a point O (fig. 1). Frank (private communication) calls 
such a meeting of dislocations a “ node’. A convenient convention, 
due to Frank, for defining the Burgers vectors of dislocations at a node is 
to look outwards from the nodal point O and describe a closed clockwise 
circuit round each line ; in the case of an imperfect dislocation the circuit 
should begin in the stacking fault associated with the dislocation. The 
same circuit repeated in the lattice as it existed before the dislocation was 
made will not be closed and the Burgers vector of the dislocation is the 
additional displacement needed to close it. With this convention the 
Burgers vectors b,, by, bs of the lines meeting at O must satisfy the relation 


3 
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If dislocation 3 is perfect it can move in any plane v where bs. v=0, 
and if neither s, nor s, lie in the plane v then s, can spiral in this plane, 
producing repeated increments of slip on it. Tt is essential, if this process 
is to produce a simple multiplication of slip on this plane, that the vectors 
b, and b, should also lie in the plane. However, relation (1) requires only 
Prat 

Vet Oia Da en a asl ke pee ee (1) 


and interesting topological consequences follow if we extend Frank and 
Read’s mechanism to the case where b, and b, do not lie in the plane v. 
There is then a component 


te-(beiyv)y=(— by. Vy) 7 well oc - (8) 


perpendicular to the plane, and the dislocation line s, is displaced by t for 
every complete revolution that it makes about s,. Hence, instead of 
repeated slip on the original plane v, a slip b, takes place on each of a 


Pig. Tt 


8 


series of planes parallel to y and spaced tapart. The sweeping dislocation 
climbs, as it were, round the helical surface that is formed because the 
Burgers vector of the pole dislocation has a component perpendicular to 
the sweeping plane. The shear deformation which results is thus of the 
homogeneous type that occurs in twinning and shear transformations. 

A mechanical twin or transformed structure generated by this process 
will build up conically in a way similar to the pyramids suggested by 
Frank (1949) in his theory of crystal growth, and the sweeping dislocation 
will rapidly assume the form of a helix, since the parts near the pole will have 
the liighest angular velocity. The formation of twins from sweeping 
dislocations connecting two poles is also to be expected. There is thus a 
close formal correspondence with the theory of crystal growth, the 
terraced steps of the latter being replaced in the present case by imperfect 
dislocations separating the transformed and untransformed material. 


2R2 
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The usual elastic interactions will occur between the sweeping dislocations 
and other singularities in the crystal, while the final shape and size of the 
twinned piece of the crystal will be determined by its strain energy and by 
the positions of large obstacles in the material. 

If mechanical twinning in a given crystal is to occur by this process, 
several conditions must be satisfied : 

1. The sweeping dislocation must produce the right shear displacement 

to generate the transformed structure on the sweeping plane. 

2. The Burgers vector of the pole dislocation must have a component 
perpendicular to the sweeping planes that is equal to the spacing of 
these planes. 

3. The pole dislocation must be anchored strongly enough to prevent it 
from moving under the stress causing the sweeping dislocation to 
move. Strong anchoring is ensured if the pole dislocation is of the 
sessile type. 

4, The sweeping and pole dislocations, together perhaps with other 
associated dislocation lines, must form a node and in this node the 
sweeping dislocation must be free to move in a sweeping plane which 
is intersected by the pole. 

In the following sections it will be shown that these conditions can be 
satisfied for {112} mechanical twinning in the body-centred cubic lattice, 
while in the case of {111} twinning in the face-centred cubic lattice the 
second part of the last condition is not satisfied, with the result that only 
monolayer twins can be formed. 


§ 3. THe Bopy-CenTRED CuBic LATTICE. 


A twinned structure can be formed in this lattice by a shear of 1//2 or 
1/2 in a (111) direction on a {112} plane, and reliable evidence exists to 
show that the Neumann bands in «-iron are mechanical twins of this type 
(e. g. Smith, Dee and Young 1928). Consider a portion of a line of perfect 


; : a P ; 
dislocation of 5 [111] type which lies in (112) but not along [111]. Then 


(fig. 2), if sufficient energy is available from external sources, a dissociation 
of the type : 


a a Cae 
SUE Sa [U2 1 TT) ee! 
can occur along a length BO of the line, forming nodes at B and O. 
Suppose for definiteness that the Burgers vectors of the lines OA, OB and 
a a 
OFEDB are, using the nodal convention of §2, 511), — 3[112] and 


7 Aaa 
= gill. Then the line OB is pure edge and sessile, while the line 


BDEO, in moving away from OB, has formed a sheet of stacking fault in 
(112). This stacking fault is, in fact, the simplest possible on the (112) 
planes if the planes are to retain the original packing of nearest neighbours, 
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and corresponds to the insertion of two fault planes in the (112) stacking 
sequence ; this sequence is a congruence to modulus 6 and the stacking 
fault is of the type 12345616123456. Since the line BDEO is glissile it 
may in its motion become pure screw ; a length OE can then move into 
either (121) or (211), since it is a twinning dislocation for these planes. 
On account of the unidirectional character of the twinning movement the 
line OE must move so that the new stacking faults generated on (121) or 
(211) make acute angles with the original fault in (112). The node at O 
can now produce twinning if OF sweeps about AOB as pole, since the 
components of the Burgers vector of OA perpendicular to (121) and (211) 


are = [121] and = [211], respectively. The sweeping dislocation thus 
a) 


Fig. 2. 


Dissociation in the body-centred cubic lattice of a unit dislocation. AO=unit 


dislocation 5 [111]; OB=sessile dislocation — 5 [112]; OFEDB=twinning 


dislocation — Fs [111]. 


climbs along OB into the stacking fault as it rotates ; it can pass through 
the latter since this is identical with a mono-layer twin on (112) and the 
twins on the three planes (112), (121) and (211) have the same direction of 
twinning movement. The movements possible at the node O have been 
described explicitly but complementary movements can also occur at B, 
again causing a macroscopic twin on (121) or (211) to grow into the stacking 
fault. 
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Or 
I 
10.2) 


§4. Toe Facr-Cenrrep Cupic Larrice. 


The analogous dissociation in this case is 
a a marie 
= —[111]+ =[112]. sbTE) ee eae oe 
5 [110] 111 ]+ 5 112] (5) 


Here, = [110] represents a unit, perfect dislocation; the line of this 


dislocation could lie in any direction and the two {111} type slip planes 
on which it could glide are (111) and (111). When the line lies in the (111) 


a 
plane the above dissociation is possible because 3 [111] represents a 


sessile dislocation in this plane. It forms the boundary of a fault, such as 
1231231312312, in the stacking of successive (111) layers on one another. 
This fault has low energy because it preserves the close-packed arrange- 
ment of nearest neighbours ; the frequent occurrence of annealing twins 
on {111} planes is recrystallized face-centred cubic metals, across which the 
stacking arrangement is of the type 12312321321, also shows that this type 
of fault is energetically preferred in this lattice. The other dislocation, 


[113 formed by the dissociation can glide in (111) and in so doing 


produces a twinning shear of the layers between which it passes. 
The strength of the component of the Burgers vector perpendicular to 


(111) is, for both the 5 [110] and 5 [111] dislocations, the same as the 


spacing a/1/3 of successive (111) layers, so that these dislocations form a 
suitable pole for twinning on (111) when this is the sweeping plane. 


a Ey 
However, the 5 [111] dislocation lies in (111) so that the sweeping = [1 2] 


dislocation must meet it as it rotates in this plane. Because of this, 
only one revolution of the sweeping dislocation can occur, even in the case 
where the undissociated part of the unit dislocation does not lie in (111) 
and the nodal point can then act as a point of emergence from the plane. 
Thus only a monolayer stacking fault can be formed. 

That repeated rotation of the sweeping dislocation is impossible may 
be seen by considering fig. 3. These diagrams show a set of (111) planes 
viewed from the side, perpendicularly to [112]. Diagram (a) shows the 


., a : : Sin ed : 
unit 5 [110] dislocation before it dissociates ; here the numbers refer to the 


stacking positions of the (111) layers and the broken vertical line is a 
boundary across which the stacking numbers are changed. By crossing 
this line from left to right the number of a layer is increased by 2 ; since 
the sequence 123123 forms a congruence to modulus 3, then 4=1 and 
5=2, etc. The dissociation occurs when the broken line is moved to the 


right, as in diagram (6). The sessile 3 [111] dislocation remains in the site 
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of the original dislocation, while the sweeping - [112] dislocation is 


situated at the lower end of the broken line. In the region between these 
dislocations the stacking sequence contains the fault 12313123. 

As the sweeping dislocation moves away it extends this stacking fault 
between the central (111) layers. By rotating about the point of emergence 
it eventually approaches the sessile dislocation again, and from the 
left-hand side. It reappears here between either the central 2 and 3, or 
1 and 2, layers, respectively, according as the undissociated section of aie 
dislocation line emerges from the (111) plane on its lower or upper side. 
Diagram (c) illustrates the first of these possibilities. The final diagram 


Fig. 3 
t ' 
2 | 2 | 
Reyes fe Es ee 
| : [111] Ln ee 
__ Se 2 eae ee Sinaia sed ep wel a <n ig 
eae senate a Ce A Tins? BO eee ee ee 
2 l 2 : ; 
es. ete all ee biel er 
3 3 BS 
Jo mera aaa Mee a 
Once ee a panmrennmnen So Eh 
(a) Davita (2) 
ie : 
3 2 2 L 3 | 5s 
Se waiah Gin 22 Sl a Tes ae aE 
SS. meena ee teers 
| 
eaMiM aio cite at te camer emmmeGenslailh <i) 2 46, 
7 ash ia - 2 ; = 
i n= eameaa a ee ee a 
ho sn ae mamma <i hed «he OE a 
ae Remecee oe te 
(c) (a) 


Stages in the formation of a stacking fault on (111) by the dissociation of 


a ei, go : : 
an 5 [110] unit dislocation in the face-centred cubic lattice. 


shows the arrangement when the sweeping dislocation has rejoined the 
sessile one, after forming a stacking fault over the entire plane. The 
sweeping dislocation cannot continue the same movement further because 
if it did, it would then be making a second circuit between the same 
layers, and would create a fault of the type 1231123 in which nearest 
neighbours are not packed correctly. 

The single rotation just described is the only possible motion for the 
sweeping dislocation. The dislocation cannot move towards the left, 
starting from its position in diagram (a), because this would create a fault 
of the forbidden type (either 2311231 or 2312231) between the central (11 1) 
half-layer and one of its neighbours. Furthermore, it cannot move in any 
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other type of {111} plane because it can only glide into a different plane 
from (111) when it is of pure screw type, 7. e. when its line lies along [112], 
and {111} planes do not intersect one another along this direction. 

It follows that this pole mechanism produces only monolayer stacking 
faults in the face-centred cubic lattice, and cannot cause the growth of 
mechanical twins. It is interesting that, while stacking faults are known 
to be formed by cold working metal crystals with this lattice, no-one so far 
has demonstrated unambiguously that mechanical twins can also be 
formed in them (Barrett 1950). 


§5. DiscussIoN OF THE PRoposED TWINNING MECHANISM. 


The question of the energy involved in the dissociation of the unit 
dislocation at the node deserves attention. In both of the dissociations 
proposed in §3 and §4, the Burgers vectors of the resultant dislocations 
are perpendicular to each other, which means that these dislocations exert 
no forces on each other from their stress fields and that the elastic energy 
of the field of the dislocations is unchanged by the dissociation. However, 
as the sweeping dislocation moves away from the other one, it trails in its 
path a layer of stacking fault, and the extra energy of this fault must cause 
a force on the sweeping dislocation. When the distance between the 
dislocations exceeds a few atomic spacings the energy of the stacking fault 
is directly proportional to its area. Thus if the energy per unit area of the 
fault is « then the force per unit length on the dislocation is also «. This 
force is equivalent to a shear stress o,=«/b, where b is the strength of the 
sweeping dislocation, acting on the sweeping plane against the direction of 
the shear displacement caused by the dislocation. Little is known of the 
magnitude of e, except that it is much smaller than typical grain boundary 
energies ; the energy of {111} boundaries of annealing twins in copper has 
recently been determined by Fullman (1950), who obtains the value 
17erg.cm.*. Taking «~1l0erg.cm.2 and b6b—10-8cm., then 
o,=10°dyne.cm.-?.. This value, which is much smaller than the 
theoretical shear strength, is not unreasonable. 

The formation of the first layer of a twin will be more difficult than the 
subsequent growth of successive layers because the addition of more 
twinned layers to an existing twin crystal does not increase appreciably the 
area of misfit between the twin and matrix. It follows that, if the stress 
needed to drive a twinning dislocation through the fields of various 
obstacles (other dislocations, faults, precipitates, etc.) in the crystal is o;, 
then the applied stress needed to start twinning is o;+0,, while that needed 
to continue the growth of a twin is only o;. Twinning resembles the yield 
phenomenon in this respect and it may be useful to define upper and lower 
twinning stresses for the start and propagation of the process. The 
characteristically rapid formation of twinned crystals, accompanied by 


discontinuities in the stress-strain relation, is almost certainly due to this 
effect. \ 
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It is well-known that twins can form extremely quickly and the work of 
Mason, McSkimin and Shockley (1948) on tin, and Millard (private com- 
munication 1950) on zinc, shows that they can form in times of the order 
of a few microseconds. It is important to consider how quickly a twin 
could form by the pole mechanism. Here we may take advantage of the 
analogy between the pole mechanism and the theory of crystal growth 
proposed by Frank (1949). When a steady spiral bas built up the process 
of growth continues by a rotation of this spiral. The spiral makes v/zl 
turns per second, where, in the present case, v refers to the velocity of the 
dislocation normal to itself and / is twice the limiting radius of curvature 
about which a section of the line can be bent by the applied stress. One 
layer is added to the twinned crystal for each rotation of the spiral. 
Following Mott and Nabarro (1948) we take /=2(u/c)b, where p is the shear 
modulus and o is the applied stress ; a reasonable value is /=10-* cm. 
The velocity v cannot exceed that of sound and we shall assume a 
limiting value v=104 cm. sec.-!. Then the time needed for the twin to 
grow by one atom layer is of the order 3 x 10-% microseconds. The time 
required to build up the steady spiral is determined by the size of the 
crystal and the applied stress. In the case of a 1 mm. crystal, and using 
the above value for J, about 1000 turns would be required. Thus the 
twinned crystal can be formed by the pole mechanism in a time of the 
order of a few microseconds. 
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SUMMARY. 


The validity of the superposition approximation used in the Born—Green 
theory of classical fluids is examined by obtaining a virial expansion of 
the pressure in powers of the density. It is found that the second and 
third virial coefficients are given correctly, but that the fourth is not. 
The discrepancy is illustrated numer:cally by considering explicitly the 
case of a gas of rigid spheres. 

The corresponding equations for a two-component mixture are also 
considered and it is shown that the equations are mutually consistent only 
when linearized. 


$1. INTRODUCTION. 


As is well known, the Born—Green theory of classical fluids embodies a 
so-called ‘‘superposition approximation’”’, first introduced by Kirkwood 
and Boggs (1942, 1949). In the notation of Born and Green, this is 
expressed by the equation 


} 


(123) ___,, (12 23 31 3 
1 1 ee eins | et ee 


where the n’s denote appropriately normalized configurational distribution 
functions. And it has been asserted (see, for instance, Rodriguez 1949) 
that this approximation is valid at sufficiently low and sufficiently high 
temperatures: 7.e. for almost perfect crystals and near-perfect gases. 
It is our present purpose to examine in some detail the validity of 
equation (1) at high temperatures or, more precisely, at low values of 
the density of the assembly. We shall, in fact, examine how far the 
Born—Green theory predicts expressions for the successive virial coefficients 
of an imperfect gas in agreement with the exact expressions for these 
coefficients that have been given by Mayer (1937-8). 
We shall find that, if the pressure is written in the form 


p= pkT(1+ Bp+Cp?+ Dp*5sc:..den)\ jp Cee ae (2) 


where p denotes the number of systems (atoms) per unit volume in the 
assembly, then B and C are given correctly by the Born—Green theory 
Se Np a el 
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but Dis not. The discrepancy will be illustrated numerically by consider- 
ing explicitly the case of a gas of rigid spheres, for which the calculations 
can be performed very simply. 

Finally, we shall consider the corresponding equations for a two- 
component mixture, and shall see that in this case a consistent formalism 
can be obtained only by subsequently linearizing the equations. 


§ 2. Born—GREEN EQUATIONS AND VIRIAL COEFFICIENTS. 


Using a slightly different notation from that of (1), we denote n,(r)/n? 
by g(r) or, more precisely, g(r, p,T). The radial distribution function, 
g(r), is then so normalized that g(r) > lasr— oo. Also, of course, n,;=p. 

Then the equation for the pressure, 7, is 


p=pkT— —~ Eat Napa le ari Or mets eee, AS) 


where ¢(7) is the interaction energy sDeesreon two systems distance 7 apart. 
This equation is exact, and does not depend on any of the subsequent 
approximations of the Born—Green theory. This theory, which makes use 
of equation (1), is essentially epitomized in a non-linear integral equation 
for g(r), namely 


log gir)+ Ger —np | | (#0) 2 g(t) —1) at ge) pp as 


(4) 


In order to obtain from (3) and (4) a virial expansion of the form (2), 
we must expand g(r, p, T) as a power series in p. Following de Boer 
(1949), we write 


2 n 
G(s Ps T)=| exp [sre | | ave-+as5 a Cae +0, 25 +... | 
(5) 


and then (4) gives a set of equations for the a,(r, T). Substitution in (3) 
then yields the requisite virial expansion. 

Writing exp [—4(r)/kT]—1=<(r), it is easily verified that the equations 
‘for a, and a, are 


a(r=—a | ec er) dia'(s) ds. . . (6) 


[ (s?—?) os a(t+r) dt a,(s)a'(s) ds 


—2r | “PP (ey a(t nlatt-+r)-+1] dt a!(a) ds. 
0 


“8 


(7) 
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Consequently, if we denote the result of substituting these expressions, 


in (3) by 
p=pkT(1-+B,p+C,p*-FDyptee ie, 
then 
12 
=2 fe! ; «| lis beanie ae 
Bi=5 |, 4rra'(r) dr, (8) 
C= 5) 4aro!(r)ax(r) DS, ate a te eg 
6Jo 
Di= pp] 477a'(r)aa(r) dr. fellatio ¢ri 
12) 4 


These formule, which give the successive virial coefficients on the 
basis of the Born—Green theory, embodying the Kirkwood—Boggs super- 
position approximation, have now to be compared with the exact formule 
for these virial coefficients which we know from Mayer’s work are 


io 6) 


B=—2n | ra(r) dr, a tie Arie Te 


0 


C=— 5] [ aretevtar dra dr, Se te ee 
and 
D=— sl{J [% 1 9% 93% 31% q1% 49% 4g 6X4 9% 93% 3 1%q1%4o 
+304 9% o3%g1%q1]ATgATgd7T,, . «. (13) 


where « ;;==«(|r;—r;|), and dr; denotes a volume element at the point r;. 
It is obvious at once, by partial integration, that B,=—B. Before 

proceeding further with the comparison, however, it is convenient to 

derive some useful mathematical formule involving Fourier transforms. 


§ 3. Some MatHEmaticaL LEMMAS. 
I. With the notation illustrated by fig. 1 


e 


| Sein 88 be Ail?) fo(7"9) bagi Sultan) fneaTnsa) dr, at, ed ai Eas dt, 


sin Ux 


3n—1 po 
=(2n) = | WF ,(u)F (tt) « . . » Foy (tt) du 
1 (oa) 
where UF (u)= Van) ip i rf (r) sin ur dr 


and the /f; are even functions of r. 


This result is easily obtained by induction. We start by considering 
the case illustrated in fig. 2. Since 


dr=r,r, dr, dr, dx/s, 
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where yx denotes an angle of rotation about AB, we have 


| falrdfalrs) dtr=27 by eS 


210 sary s 


dr, dr.. 
1 2 : 
Writing 7,f,(7,)= V(@n) 7 uF ,(u) sinr,judu, and interchanging the 
‘77 —o 


Fig. 1. 


A 


order of integration, we find 


To 


n su ; 
[Al 11)fel’2) d7= SOO aa |. d a UF {(U)fofo(7o) Sin Tou du 


gin su 


du. 


=o F uF ,(u)F o(u) 


Considering next the integral (see fig. 3) 


T= || flrfalradfalrs) dr de’ 
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we have 
rifi(7s)rofel”2) 8173f3(1'3) ds dry 
t 


T=(2n)?{.... | dry dr, . 


va co rt+1s i 
LOA | | | uF ,(u)F (uw) — rafalt'y) ds dr'y du 
—a/ (0 t—?, 


sin tu 
t 


du. 


=(2ns® | uF (uF (uF ou) 
The general result stated in I. follows by simply repeating the argument. 
Letting «— 0, we have 


II. With the notation illustrated by fig. 4 


[. tail: Jula)falrs) oe Snei(Tnt1) Ity It,» -- dT, 
, 3n—1 7a 
=(27) 27 | WF (UyF (u) .... Fyy4(u) du. 
Fig. 3 
ake 
ue 


ats 


ae | 
we 


£ 


These results are also proved in a paper by Montroll and Mayer (1941), 
but their method of derivation is considerably longer than ours. For the 
formule to be valid the interatomic potential must decrease more rapidly 
than for purely ionic interactions, which is certainly so for the van der 
Waal’s forces normally responsible for gaseous imperfections. 


§ 4. Discussion or C anp C,. 
On account of II. above, equation (12) gives at once 


O=—}(2n)52 is w[B(u)P® du, . . . . . (14) 
where ' ; 
up(u)= heise iy ") sin wr dr 
ae pres) sin ay ee Parte ee rl CES Y) 


a(7r) being regarded as an even function of r. 
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To obtain from (9) an analogous equation for C, we must first consider 
a,(r). On replacing (¢+-r)x(t+7) in (6) by its Fourier transform according 
to (15), namely 


i co 
a) (2m) [ p uB(w) sin u(t+r) du 


a integrating first with respect to ¢ and then with respect to s, we find 
that : 


ray(r)=2n | ULNA SUNG OU. 0. eo ee LG) 
Consequently, by (9), 


, 4n? 
C= 5 rx(r) | u[B(u)]? sin ur du dr 
< 0 — 
Fig. 4. 
dt 


Integrating partially with respect to r we have, since «(r) vanishes at 
infinity, 


qe re 
=f” {veneers + YS” waza g,(wBle)} de 


and the second term on the right vanishes since for interaction potentials 
of the type under consideration uf(u) vanishes at infinity. Thus 


C,=—H2n)s2 [wrt B(w)]? du=C 


and the third virial coefficient is given correctly by the non-linear integral 
equation of the Born—Green theory. 
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§ 5. Discussion oF D anp Dj. 


On account of IT., and by applying Parseval’s theorem to equation (16), 
we have for the last term on the right hand side of (13), 


mt ie ih a) 
leer lH] 12% og%gq%q1 IT, ATs dry=—42n)*| u*[B(u)|* du 
=— S| air) dr. a 
0 


For the second term on the right hand side of (13), application of I. and 
II. leads to 


Ff f ararcetastartss dri dre drp=—2m5 [wR (w)[Blu)} de, 


— © 


where 


UF (u)= onl. if(t) sin ut dt and tf(t)=27«(t) pe ul B(w)]* sin ut du. 


‘Using Parseval’s theorem, and (16), we thus obtain 


- iI | 11 o%o9%q1%41%49 ATy AT, dt4= —3rr IE rea(rjaz(r)dr. . . (18) 
We have not been able to obtain any such simple formula for the first 
term on the right hand side of (13), but we shall find that the expressions 
involving a@,(7) in (17) and (18) arise naturally in the evaluation of D,, 
to which we now turn. 
From equations (7) and (10) we obtain, after partial integrations of (7), 
7 ive) 


Die - | b ra’ (r)az(r) dr (i.) 


pe a 0 72a" (r) (hn [ (s?@—#?)(t-+-r)a(t-+1r) dt . a,(s)«(s) ds dr (ii.) 
0 0 = 


oS ra'(r) i [ena dt . sa,(s)a(s) ds dr (iii.) | (19) 
qq 5 ears 

sbi “acon | | (br jalt+njay(t++n) de. x(a) ds dr (iv. 
Aq? (% oO rs : 

ne I. ra (r) I, | a (¢-+r)a,(t-+1r) dt . sa(s) ds dr. (v.) 


It is convenient to number the successive terms on the right hand side of 
(19), (i.), (i), ... (v.). 
Then using the techniques already illustrated above, we can show 
without difficulty that 
ony 


(v.)=— Oh raz(r) dr=(17) 


~*~ 


ee eee oe - Se ee 


EEE eee 
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and 
= ra(7)a2(r) dr+ ie ra(r)a,(r)a,(r) dy 
J0 


Gi=ativ == | 
so that 


0 


(i.) + (iii.) + (iv.) = — Br | rax(r)a2(r) dr—=(18). 
0 

Thus the second and third terms on the right hand side of (13) are 
accounted for by the terms (i.), (iii.), (iv.) and (v.) of (19). But the remaining 
part of D,, namely term (ii.) in (19), is not equal to the remaining part of 
D, i. e. to the first term on the right hand side of (13). This is clear from 
dimensional arguments and, as we shall see, is confirmed when we make 
explicit calculations for the case of hard spheres in §7 below. For the 
purposes of calculation it is convenient to write (ii.) in the form 


: Qn® po dv .df1 [pup . 
ae ee 2! eee (; J | i URE ars) f on D) 
(i1.) 3 i? f ay (r)a(r) = E ml yn ur ant dr. (20) 
We thus see that the Born—Green non-linear integral equation, (4), 
for g(r), deriving from the superposition approximation, fails to give 
correctly the fourth virial coefficient, D, of an imperfect gas; and that 
the reason for this failure is that it does not do justice to the term involving 


JJ §¢40%29%34%41%q0%43 AT. dT, a7, in the correct expression for this virial 
coefficient. 


Before giving the results of explicit calculations for the simple case of 
a gas of hard spheres it will be well to consider how these results are 
affected, if at all, if we base discussion not on (4) but on a linear approxi- 
mation to this non-linear equation. Most of the attempts made so far 
to determine explicit distribution functions for a given intermolecular 
potential have been based on such a linearization technique. We confine 
our discussion to the particular method of linearization used by 
Rodriguez (1949). 


§ 6. Errect oF LINEARIZING THE EQUATION FOR g(7). 
Rodriguez has shown that if, in (4), we write 
g(r)=exp [f(r)—¢(/AT], . . 2... (21) 


keep only terms linear in f(r) and replace certain products f(r)a(r) and 
f(r)a'(r) by (e—1)a(r) and (e—1)«‘(r), respectively, where 


= Me oi des egal 22) 
| ra(r)dr 
Jo 


then the resulting linear equation can be solved, the solution being given by 


1 %° ¢s[B(s)]*? sin rs 
Ory Ne = eee 


SER. 7, VOL. 42, NO. 329.—JUNE I9Q5I 25 


(23) 
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where A~!=(27)39. The precise nature of these approximations is not 
at all clear, but our final paragraph will suggest that when dealing with 
binary mixtures some such technique is not only convenient but necessary, 
and it will be well therefore to discuss briefly the virial coefficient problem 
on the basis of equations (23), (22), (21) and (3). 

The procedure is exactly as before, except that now we must also 
expand « as 

e=1-+e,pcgp?4 s5- > 2) 2 | gl tee eee 

There is no need to enter into the details of the calculation. It is sufficient 
to state that if we express the pressure p as 


p=pkT[1+Byp+Cop?+Dop?+ « - - J. 


Then we find 
Beeb — 5, 
G.=G=c, 
but 
8 2 ioe) oa) 
Di > Cy | 72a’ (7) { s[8(s)]? sin rs ds dr (a) 
0 —o 
2a ng 2 fis Ms 
- ; (2m)5? | ra'(r) | s[8(s)]}? sin rs ds dr (b) | (25) 
0 —o 
4773 hoe) foe) 2 
+ =r ra’ (r) | { s[B(s)]? sin rs as dr, (c) 
0 — 0 
where 


C;=(27)82 |. ut B(u)]> du] | 1%a(0 dr=30/2B. i sues 
J 0 


It is then not difficult to show that, comparing the terms of (25) and (19), 
(b)=(v.), 


bat (a) A(ii.) + (iii.) + (iv.). 
In fact (a)=3C?/B. «oe as es et 

Thus linearization does not affect the second and third virial coefficients 
B and C, but does affect the fourth virial coefficient. 

We shall now illustrate these results by using the above formule to 
calculate B, C, D, and D, for a gas of hard spheres. The value of D (as 
also B and C) for such a hypothetical gas has already been calculated by 
Happel and by Majumdar: see Fowler and Guggenheim (1939). 


§ 7. CALCULATIONS FOR A GAS OF HaRD SPHERES. 
We consider hard spheres of diameter a. Then ¢(r)= 00, r<a and 
¢(r)=0, r>a. Whence 
a(r)=—1, O<r<a; a(r)=0, a<r. 
‘Then 
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and, by (16), 
ra,(r)=0, |r|>2a 


7 
Se r(r—2a)*(r+4a), |r|<2a. 
From (11), 
27a 
B= 3 =4v=b say, 

where v is the volume of one of the spheres. 

From (14), 
5728 
— —5 2 
is 3b?, 


We now consider the various contributions to D, and D,. They are 
given in Table I. below : 


TABLE I. 


D, Value D; 


(i.)-+(iii.)+-(iv.) | 6847 vans 
from (18) 4480 b3 (a) from (27) 


35 


(v.) from (17) — 345: (b) from (17) 


(ii.) from (20) = ean b3 (c) from (i.) 


Total D, ed bs Total D, i 


= 0:3958 68 


These values have to be compared with the value D=0-286963 calculated 
by Happel and by Majumdar. We see that the non-linear integral 
equation leads to rather too small a value for D, while Rodriguez’s 
method of linearization yields too large a value. 

Finally we would observe that (26) shows that for hard spheres « should 
increase with the density p. If, however, we adopt the crudest form of 
linearization and take «=1 throughout, then B, and C, are unchanged, 


but D, loses the contribution (a) and becomes Dj=—#2 3, which 
has the wrong sign. Here we have not ignored f? in finding p; if we do 
so we get Dj=—# b3, which is also negative. 


§ 8. Brnary MIXTURES. 


In this discussion we shall try to keep the mathematical exposition as 
short as possible since the calculations though essentially straightforward 
are necessarily cumbersome if presented at all fully. To obtain the 


282 


592 G. S. Rushbrooke and H. I. Scoins on Virial 


equations we have simply to follow Born and Green’s procedure. We 
derive first an infinite set of equations for the equilibrium distribution 
functions, and then cut these off by means of the superposition approxi- 
mation. We now need, however, four equations of type (1), namely 
(with a self-explanatory notation) 


My ya(123)=My 4(12)M4.4(23)r4 4(31)/M4, (28) 
1 4p(123)=My 4 (12)r4 p(23)4 p(13)/ning, 
and the two equations which result from interchanging A and B. 


We then obtain, in place of (4), four non-linear integral equations which 
may be written 


ia (80) galt) at gyal SS ds 


roo 78 t 
trpn| [ (—t)*fginlt+r)—1] dt gyn) 
. ds, Oren 20) 


where i=A or B, j=A or B, and the order of the suffices is immaterial. 
In terms of the n’s, g;;(7)=7,,(r)/nn; and > 1 as roo. Moreover 


| my=N,/V=p, and npg=N;/V=pz. 
The equation for the pressure is 


p / 
Ev ae fl [eg a(r)b44(7) +202 p94 8 (7)h4p(7) 


+x pn(")Opp(")]}4rr3 dr, . . . (30) 
where p=(N,+Np)/V, (s=N,/(N4+Np,) and 4, =Npz/(N,+Np). 
To obtain a virial expansion we proceed exactly as before except that 


now we have to use double power series in the two variables py and pz 
instead of simply power series in p. If we write 


p=pkT[1+Bp+Cp?+Dp?+ ...] 


then we find 
B=—27 E ib aa a(r)r? dr+-27 0p 4 p(r)r? dr+ay ip opp)? ar] 
“0 0 
31 
and Pi 


rf ; 
C=— 3 E i 4 4 (712)%4 4 (%o3)%4 a (713) drs drs 
+3052, {| 4 (11 2)%4 B(7o8)%4 B(I'13) dr, drs 
+32 ,x% {J Ap Bll12)%4 B("93)%4 B("1g) AT. ATs 


+28 {J @3p(P12)%pp(es)tp p(s) dT» ars | een C29) 
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And these are the correct expressions for the second and third virial 
coefficients of a binary mixture (Mayer 1939). 

To this stage the calculations, though laborious, are essentially straight- 
forward. The only mathematical tool required, apart from algebraic 
manipulation, is the result 


I, is (s°—#)(t-+-r) f(t+r)g'(s) dt ds 


olteul (s?@—#)(t+r)g(t+r)f'(s)dtds, . . . (33) 


which is easily proved by introducing the Fourier transforms of rf(r) and 
rg(r). We need (33) in order to show that the two equations of type (29) 
for g,,(7), obtained by taking i=A, j=B and i=B, j=A respectively, 
are mutually consistent as far as terms linear in p are concerned. 

When however, we proceed to calculate D we meet an entirely new 
difficulty. For we now find that these two equations for g, ,(7) are not 
mutually consistent when we include terms quadratic in p. Indeed, it 
can be shown that the reason for the inconsistency can be traced to 
precisely the same defect in the superposition principle that led to D,~¢D 
for a one-component assembly. 

We thus find that when the Born-Green formalism is extended to binary 
mixtures the superposition «pproximation leads to four non-linear 
equations for the three equilibrium radial distribution functions, and that 
these four equations are not mutually consistent. If, however, we 
linearize the equations in the manner due to Rodriguez, then in place of 
(29) we have the four equations 


rfig(?)=27p 4 | 


‘ 


[o) 


[CEP eg +r) Heraaialt$r)] Ut e548%54(8) ds 


-o r 


| 2 ((-+7)[fiplt +r) +€:3%p)t+7)] dt ¢; p80; ,(8) ds, (34) 


0 


i=-A, B and j=A, B; and it is easy to show that these equations are 
entirely mutually consistent. Their solution, by the method of Fourier 
transforms, can be found without difficulty : we shall not, however, give 
it here as the expressions are rather unwieldy and we do not wish to 
discuss them further at present. 


REFERENCES. 


pE Bour, J., 1949, Reports on Progress in Physics, 12, (London: The Physical 
Society), p. 305. * 

Born, M., and Green, H. 8., 1949, A General Kinetic Theory of Liquids, 
(Cambridge : University Press), Ch. I. and II. ; 

Fowter, R. H., and Guecenuerm, E. A., 1939, Statistical Thermodynamics, 
(Cambridge : University Press), § 717. 

K irxwoop, J. G., and Boaas, E. M., 1942, J. Chem. Phys., 10, 394. ' 

Maver, J. E. et al., 1937-8, J. Chem. Phys., 5, 67,74; 6,87; 1989, Tbid., 48, 71. 

Monrrott, E. W., and Mayer, J. E., 1941, J. Chem. Phys., 9, 626. 

Ropriausz, A. E., 1949, Proc. Roy. Soc. A, 196, 73. 


[ 594 ] 


LXIL. The Critical Magnetic Fields of Aluminium, Cadmium, Gallium 
and Zinc. 


By B. B. Goopman and E. MENDOZA *, 
Royal Society Mond Laboratory, Cambridge 7. 


[Received March 7, 1951.] 


ABSTRACT. 


The low temperatures produced by the adiabatic demagnetization of 
a paramagnetic salt have been used to cool specimens of aluminium, 
cadmium, gallium and zine, and their critical fields have been measured 
down to 0-1°K. The magnetic measurements on the metal and the 
salt were made independently and the experimental arrangement had 
the further advantage that the quality of the thermal contact between 
the metal and the salt could be studied. The results could be accurately 
expressed by relations of the form H,=H,(1—(T/T,)?); values of Ho 
and T, for each metal are tabulated. 


§ 1. INTRODUCTION. 

CoMPARATIVELY little information is at present available on the critical 
fields of superconductors with low transition temperatures. In this 
paper measurements of the critical fields of aluminium, cadmium, gallium 
and zine are described. The low temperatures were produced by the 
adiabatic demagnetization of a paramagnetic salt, using a technique 
already described briefly by one of us (Mendoza 1948), in which the 
metal specimen and the salt were separated but were in thermal contact 
through a long copper rod. This arrangement has the advantage that 
it is possible to make magnetic observations on the metal specimen and 
the salt independently. Moreover the condition of the specimen can be 
more easily controlled than in the method used by Kurti and Simon 
(1935 b) and more recently by Daunt and Heer (1949 a, b), where the 
salt and chips of the metal were pressed together into a pill. The 
superconducting transitions in various steady magnetic fields were studied 
by measuring the effective susceptibility of the metal in a low-frequency 
alternating field ; the specific resistance of the normal metal could also 
be deduced from these measurements. 


§2. Tur ExpErRtMeNTAL ARRANGEMENT. 

A diagram of the low temperature part of the apparatus is shown in 
fig. 1. The salt-tube T stood in a Dewar flask of liquid helium, itself 
surrounded by a Dewar of liquid air. The metal specimen S was gripped 
in the cup at the lower end of the copper rod C. Strips of copper foil F 
with a total area of about 6 cm.? were hard soldered to the top of C: 


* Now at the H. H. Wills Physical Laboratory, Bristol. 
+ Communicated by Dr. D. Shoenberg. 
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there was no evidence that the magnetic observations were affected by 
the hard solder. A paste of powdered paramagnetic salt (either ferric 
ammonium alum or chromium potassium alum) mixed with a solution 
of a plastic cement in acetone was pressed around the foil in a cylindrical 


Fig. 1. 
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The low temperature part of the apparatus. 


die under a pressure of 2000 atm.; the resulting pill K was 1-6 cm. in 
diameter and about 4cm. long. It was found that this arrangement 
gave adequate thermal contact between the salt and the superconductor, 
and that within five minutes of demagnetizing the salt to 0-1° K. good 
thermal equilibrium had been reached. The amount of heat reaching 
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S and C from stray sources was reduced by a cylindrical shield of copper 
foil H in good thermal contact with the salt K. Vertical slots were cut 
in H and in the cup holding S in order to reduce the effect of eddy currents ; 
the slots in the cup also made it easier to grip S. The whole system was 
suspended by cotton threads L from the framework G which stood inside 
the salt-tube T. Hot gas molecules coming down the pumping tube P 
were prevented from impinging on the cold salt K by the radiation 
shutter R. 

The longitudinal magnetic field needed for the critical field determina- 
tions was provided by a solenoid which could be suspended around the 
liquid air Dewar. This solenoid gave a field which was uniform 
within } per cent over a length of 3-0 cm., and the effect of its field on the 
temperature of the salt was always negligible. 

M, and M, are mutual inductances surrounding the salt and the 7 
superconductor respectively. Since the value of the empty mutual 
inductance was not known very accurately, only changes in susceptibility 
of K or S could be deduced from measurements of the mutual inductance. j 
Owing partly to losses in the circuit and partly to losses in the specimen ; 
(particularly for the metal in the normal state) each mutual inductance 
can be regarded as complex and will be denoted by M’—jM”"’. Values of 
M’ and M” for both mutual inductances were measured by a Hartshorn 
A.C. bridge. The measuring current was supplied by a 40 c./s. frequency- 
stabilized alternator and detection was by means of an amplifier and a 
vibration galvanometer. 

The contribution to M; was of order (400/T) wH from the salt and the 
change of Mj corresponding to the transition of the metal from the normal 
to the superconducting state was of order 40 .H. By using a measuring 
field of amplitude | gauss the bridge could be set to an accuracy of 
0-05 uH in the absence of a steady magnetic field, but only to 0-3 wH ina 
field owing to the effect on the galvanometer of fluctuations in the field. 


lg ay 


TABLE I, 
Characteristics of the Metal Specimens. 


Estimated Length Diam. 


Boures % purity (mm.) (mm.) (p/P ars) X 10° 
Aluminium  J1011 99-995 26-2 6-5 0-074 18 
Cadmium N 99-996 28-5 71 0-074 3-1 
Gallium J3419 ~99-99 26-8 3:00 = 0-023 1-0 
Zine J 1790 99-999 27°5 6-8 0-074 14 
J=Johnson, Matthey & Co. N=New Jersey Zinc Co. 


The quantities n and p/p.73 are referred to later in the paper. 


§ 3. THE SPECIMENS. 
The characteristics of the metal specimens are given in Table I. The 
aluminium and zine specimens were both turned to an accurately 
ellipsoidal shape. The cadmium and gallium specimens were cast in 
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thin-walled glass tubes which were approximately ellipsoidal in internal 
shape, and the glass was afterwards cracked away. Each specimen was 
annealed in a helium atmosphere to remove internal strains, and then 
lightly etched in acid to remove surface impurities. 


§ 4. EXPERIMENTAL PROCEDURE. 
(a) The temperature of the specimen. 


The demagnetizations were made from magnetic fields of up to 
20 kilogauss provided by a water-cooled solenoid which could be raised 
into position around the liquid air Dewar. It was found that unless 
demagnetizations were carried out from a helium bath temperature of 
about 1:5° K. or more the rate of warming of the cold salt was rapid 
above about 0-7° K. It seems probable that at lower initial temperatures 
the out-gassing of the salt tube was much less complete, and after 
demagnetization helium from the salt-tube wall condensed on the cold 
salt, only to evaporate and destroy the thermal vacuum when the salt 
warmed to about 0-7° K. Under the best conditions the natural heat 
input into the cold system was about 3 ergs sec.~}. 

The temperature of the salt after demagnetization was determined by 
standard procedure (Kurti and Simon 1935 b). The mutual inductance 
of M, was first measured at three or four temperatures in the helium 
range, while helium gas in the salt-tube at a pressure of about 0-1 mm. 
of mercury maintained the salt in thermal equilibrium with the helium 
bath. The helium bath temperature was deduced from the pressure in 
the cryostat using the 1949 scale of temperature (van Dijk and 
Shoenberg 1949). A graph was plotted of M; against 1/(T—4), where 
A=f(4r7/3—N); f is the filling factor, 4 the Curie constant per ¢.c. and N is 
the demagnetizing factor of the external shape of the pill ; the parameter 4 
which was usually about 0-02° K. was introduced in order to allow for 
the non-spherical shape of the specimen. The graph was a straight line 
as was to be expected for a salt obeying a Curie law. After demagnetiza- 
tion the value of M’ was again measured and a magnetic temperature T* 
was deduced from the value of 1/(T—A) given by linear extrapolation 
of the graph. The values of T* obtained at low temperatures differed 
slightly from thermodynamic temperatures T owing to departures from 
Curie’s law and were corrected on the basis of various published data 
(Casimir, de Haas and de Klerk 1939, Cooke 1949). The reliability of 
this procedure for deducing T is indicated by the agreement between 
estimates of T for a given critical field of a zinc specimen using ferric 
ammonium alum and chromium potassium alum respectively. At 
0-2° K. the two sets of measurements were in agreement to better than 
10-2 deg. although the values of T*—T were of order 3x 10 deg. 


(b) Measurement of the critical field. 


After the demagnetization the specimen was superconducting and 
measurements of Mj and Mj’ were made as the magnetic field of the 
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solenoid B was steadily increased. The real part M3 was constant as 
long as the specimen was wholly superconducting, but owing to coupling 
with the solenoid there was a small increase in Mj’ as the resistance of 
the solenoid circuit was lowered. At the entry to the intermediate state 
a sharp rise of both Mj and Mj’ was observed, and the steady field 
corresponding to this was assumed to be (1—m) times the critical field H, 
of the superconductor, where 47n was the demagnetizing coefficient of 
the specimen. The cadmium and gallium specimens were not accurately 
ellipsoidal but the resulting uncertainty in allowing for the demagnetizing 
effect was only of order } per cent. Immediately after the sharp rise 
in Mj and M3’ had been observed the value of Mj; was measured, thus 
giving the temperature at which H, had been measured. Owing to the 
gradual warming of the salt a small correction was necessary to allow 
for the rise in temperature between the measurement of H, and the 
measurement of M/. By repeating the whole procedure as the salt 
warmed up to bath temperature a series of values of H, at various 
temperatures between 0-1° K. and the transition temperature was 
obtained. 

The transitions in the absence of an applied magnetic field were 
studied by following M3 and M3’ as the metal warmed through the 
transition temperature. Values of T, deduced in this way agreed well 
with those obtained by extrapolation of the critical field curves. From 
the change of Mj and M3’ between the superconducting and normal 
state the value of the specific resistance of the normal metal could be 
calculated using standard eddy current theory. For this purpose the 
specimens were treated as if they were long cylinders ; the error due to 
this assumption is probably not more than a few per cent. Values of 
the ratio p/ps73 of this resistivity to the resistivity at 273° K. which are 
a guide to the quality of the specimens, are given in Table I. 


§ 5. CORRECTIONS AND SOURCES OF ERROR. 


The measurements of the specimen temperature and the corresponding 
critical field were subject to a number of corrections and sources of 
error which need discussion. 

(a) Temperature. 

In preliminary measurements by one of us (Mendoza 1948) it was 
found that the thermal contact between the copper rod and the salt was 
inadequate for the amount of heat being developed in the rod by eddy 
currents (no heat of course was developed in the superconductor). This 
was revealed by the fact that the values of T recorded for a given H, 
inc eased when |a ge measuring currents were used. In the present 
arrangement the thermal contact between the copper rod and the salt 
was improved ; a correction was applied for the temperature difference 
between the rod and the salt produced by eddy current heating but this 
was only of order 10-? deg. for the measuring currents actually used. 
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There remains the possibility that the temperatures of the specimen and 
the salt may differ owing to a thermal gradient set up by the natural 
warming of the whole system. Estimates of the probable value of the 
temperature difference between the salt and the specimen due to this 
cause showed that under the actual experimental conditions and at 
temperatures above 0-1° K. it was unlikely to exceed 10-3 deg. and no 
correction was therefore applied. 

Taking into account these possible sources of error and possible 
inaccuracies in the 1949 scale of temperature, the final values of the 
temperature of the specimens can probably be considered reliable to 
within 10-? deg. 


(6) Critical yield. 

The effect on the earth’s field was allowed for by making the 
measurements with the solenoid field alternately with and against the 
vertical component of the earth’s field. In order to allow for the finite 
amplitude of the measuring field the peak value was added as a correction 
to the applied field. Measuring fields of between 0-5 and 5 gauss r.m.s. 
were used in the experiments and after the application of this correction 
results were obtained with an internal consistency of better than 0-5 gauss. 

The main uncertainty of interpreting the critical field measurements 
lies in assessing the effect of strains and impurities in the specimen. Since 
the value of M, was determined by the currents flowing in the specimen 
and not by its magnetic moment it follows that the present method 
studies the electrical resistance of the specimen rather than its 
magnetization, and the effect of strains or impurities in the specimens 
is to give values of H, which are larger than the true values. An estimate 
of the magnitude of this effect is provided by the behaviour of the zinc 
specimen. The effect of annealing was to reduce H, at 0-2° K. by only 
4 per cent while the breadth of the zero field transition was reduced 
from about 410-2 deg. to not more than 5x10-3deg. If a single 
mechanism is responsible in a superconductor for the broadening of the 
zero field transition and for raising the critical field above its true value 
then it seems probable that systematic errors in the values of H, for the 
annealed specimen of zinc due to this uncertainty do not amount to 
more than about 2 per cent. All the other specimens had zero field 
transitions at least as sharp as the zinc specimen, and it seems plausible 
to assume that errors in H, values for those metals due to strains and 
impurities are of the same order of magnitude as for zinc. 


§ 6. RESULTS. 


It was found that the temperature variation of the critical fields could 
be represented by the formula H,=H)(1—(T/T,)”) ; this is illustrated 
in figs. 2 and 3 where H, has been plotted against T”, and it can be seen 
that the points lie well on straight lines. Values of Hy and T, deduced 
from these graphs are given in Table II. 
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Measurements on aluminium and gallium have previously been made 
by Shoenberg (1940) and appear to differ appreciably from the present 
results. The discrepancies are much reduced however when due account 
is taken of the differences between the 1932 and 1949 scales of temperature, 
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and of the fact that no allowance was made in Shoenberg’s measurements 
for the pressure drop between the liquid helium level and the top of the 
Dewar flask*. In fact, the diameter of his flask was only 9mm. and a 
pressure drop of order 0:07 mm. of mercury could have been expected 
at the lowest temperatures, corresponding to an error in temperature 
determination of order 3x10-?deg. Thus Shoenberg’s measurements 
may be regarded as supporting the present results. Kurti and Simon 
(1935 b) obtained in measurements on cadmium T.=0:54+0-02° K. and 
an initial slope of the critical field curve of order 100 gauss deg.-!; the 
present values of T,—0-560° K. and 103 gauss deg.—! respectively are in 
good agreement. 


TABLE IT. 
Summary of the Results. 
Te Hy “yx 104 = x 108 
° asl =) c 
rr.) (gauss) cal. mole.—! deg. Fale SIS 
(a) (0) (a) 
Aluminium 1-197 106-0 2-95 3:48* 24-9 
Cadmium 0-560 28-8 1-28 = 17-9 
Gallium 1-103 50-3 0-91 — 7-2 
Zine = 0-905 52-5 1-16 1-257 14-1 
L-5t 


y is the coefficient of T in the expression for the specific heat of the normal 
metal, and V is the atomic volume. 


(a) Calculated from the present results. 

(b) Calorimetric measuresment : * Kok and Keesom (1937); + Keesom and 
van den Ende (1932) recalculated by Silvidi and Daunt (1950) ; = Silvidi 
and Daunt (1950). 


The values of H, obtained by Daunt and Heer (1949 b) using Kurti 
and Simon’s method are of order 4 gauss lower than the present results 
for aluminium and are of order 8 gauss higher than the present results for 
zinc. It therefore seems improbable that a single factor is responsible 
for these discrepancies. It is possible that the stresses the metal is 
subjected to when pressed into a pill with the salt may be responsible 
for the high values of H, recorded for zinc. Furthermore, the difficulty 
of interpreting the complex magnetic behaviour of the mixed pills of 
salt and metal used may possibly be responsible for unsuspected errors 
in measuring both H, and T. 

The critical field measurements may also be compared with various 
calorimetric data. The difference in entropy per unit volume between 
the normal and superconducting states is given by the thermodynamical 


relation 
H, /dH, 
oe aoa \ ig ie 


* This possibility was pointed out to us by Dr. Shoenberg. 
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The experiments of Keesom and Kok (1934) on thallium and of Keesom 
and van Laer (1938) on tin suggest there is no linear term in the specific 
heat of a metal in the superconducting state. Therefore for those metals 
for which a relation of the type H,=H,(1-++(T/T,)?) holds the linear 
term in the specific heat of the normal metal must be of the form 
(VH2/27T?)T erg mole! deg.!, where V is the molecular volume. 
Values of y, the coefficient of T in this expression, calculated from the 
present results are in fair agreement with calorimetric data as can be 
seen from Table II. 

Evidence for an empirical correlation among superconductors has 
recently been put forward by Daunt (1950). He observes from a plot of 
y/V and T, that the elements for which reliable information is available 
fall on one of two straight lines through the origin. In other words he 
suggests that y/V=KT, where K=30X10~-® cal. cm. deg.-* for the 
hard group of superconductors and K=6-7 x 10~® cal. em.~* deg.-? for 
the soft superconductors. The values of y/VT, in Table II. calculated 
from the measurements on the four soft superconductors studied in this 
paper suggest that the correlation is at best only approximate. 
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SUMMARY. 


This instrument serves to obtain the size distribution of electric pulses 
coming, for example, from an ion chamber or counter. Each pulse causes 
a small steel ball to be propelled along an inclined board; the ball 
describes a parabolic path and lands in one of 30 parallel grooves. 
As balls accumulate in the grooves a histogram of the pulse size 
distribution is built up. 

The paper describes the mechanical construction and the associated 
electronic circuits. The latter serve, among other things, to eliminate 
pulses which are either too large or too small to be recorded, or which 
follow too close upon the last recorded pulse. 


INTRODUCTION. 


A KICK-SORTER (or pulse size analyser, also pulse spectrograph) serves to 
obtain the amplitude distribution of electric pulses, e.g. from an ion 
chamber or porportional counter. The pulses must first be amplified by 
a standard pulse amplifier (see Sands and Elmore 1949) ; the kick-sorter 
then sorts them into a number of channels (size groups), recording the 
number in each channel separately. 

It occurred to one of us (O. R. F.) in 1947 that a mechanical system 
might be cheaper and more reliable than the usual electronic gear which 
involves a great many valves (usually several per channel). In the 
meantime electronic kick-sorters have been much improved (see in 
particular Wilkinson 1950, Hutchinson and Scarrott 1951), but the 
instrument described here (in use since 1949) has still some advantages, 
chiefly cheapness and the fact that the information it contains is displayed 
as a histogram all the time. Its main disadvantage is low storage capacity 
(about 100 per channel) and slowness: pulses are sorted at most at the 
rate of about 6 a second. However, in many cases that is enough. 


GENERAL DESCRIPTION. 


Each pulse suitable for recording (see later) causes a steel ball (} in. 
diameter) to be propelled along an inclined board (1) with 30 parallel 
grooves (fig. 1). The ball describes a parabolic path and lands in one 
groove, depending on its initial speed which in turn depends on the 


ee eee eee ee ee ee ee ee 
* Communicated by the Authors. 
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amplitude of the pulse. As balls accumulate in the grooves a histogram 
of the pulse size distribution is built. Immediately after one ball has 
been fired a new one is deposited by an electrically operated “ feeder ” (2). 
The balls are projected by a little “ billiard cue ” (3) mounted on a coil 
which is driven like the coil in a dynamic loudspeaker. The current 
pulse which makes the coil move, and the one which actuates the feeder, 
come from the electronic pre-sorter (4) which contains some twenty 
valves. Here the pulses issuing from the standard amplifier are first 
lengthened (4.1) to 100psec.; if they are bigger than an adjustable 


Fig. 1. 


bias B their excess over B is then further amplified (4.2). The next 
section (4.3) selects pulses whose amplitudes lie within a certain range ; 
the limits of this range can be adjusted so that they correspond to the 
first and the last groove on the board. The waiting-time unit (4.4) 
closes a gate after each pulse, for a time T (0-1 to 0-2 seconds) which is 
adjusted to make sure that the next pulse does not arrive before the next 
ball is in position to be fired. The output unit (4.5) lengthens the pulse 
further and gives it enough power to drive the billiard cue; it also 
produces a suitable pulse to actuate the ‘‘ feeder ”’. 
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DETAILED DESCRIPTION. 


1. The Board (see figs. 1 and 2) was made of polystyrene, 23 cm. by 
48cm., jin. thick. Thirty grooves just over }in. wide and in. deep 
were cut in with a miller; this tapers the upper end of the groove and 
allows the ball to enter without bouncing. The board is mounted on 
three adjustable legs and is tilted by about 10 degrees in such a way 
that the line AB remains horizontal. The balls are propelled parallel 
to AC and hence describe parabolas each of which enters the appropriate 


Vig. 2. 
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groove without a break. If a ball describes a parabola which leads to 
the ridge between two grooves, then it will be caught in the lower groove 
a little lower down. 

At the lower end the grooves are closed by a removable metal bar. 
As more and more balls accumulate in the grooves, they build up a 
histogram of the pulse size distribution. Lines corresponding to 5, 10, 
15 etc. ball diameters are drawn on the board so that the number of balls 
in each groove can be quickly read off after a run. This done, the metal 
- bar is removed, and the balls run out into a gutter and from there into 
a beaker, ready to be used again. 


SER. 7, VOL. 42, NO. 329.—JUNE IQ51 2 


| 


606 S. G. F. Frank, O. R. Frisch and G. G. Scarrott on a 


On dry days the top ball in a groove often comes to rest several milli- 
metres above the remainder, presumably due to electric charge acquired 
by friction. Such a charge may slightly influence the path of the ball 
though we have no evidence that it does. A board made of metal rather 
than polystyrene would eliminate this. 

2. The feeder is a simple escape mechanism, connected to the armature 
of a small electromagnet (a commercial relay). When the armature is 
attracted the lowest ball is released; on the return motion the whole 
line of balls in the feeding tube is allowed to advance one unit, making 
a fresh ball ready for release. On being released, the ball drops by about 
14 diameters, on the spot from which its predecessor has just been fired. 

To prevent bouncing, the platform on which the ball drops is a mica 
membrane cemented over a metal ring and smeared thickly with grease 
from below. This has proved very durable and is remarkably bounce- 
free ; balls are stopped dead even if dropped from a height of several 
inches. There is some slight bouncing from the two locating posts 
which prevent the balls from rolling off. This could probably be 
improved by suitable design of the posts, and then a waiting time of 
0:05 seconds between balls should suffice. As it is, we find operation 
getting erratic if the waiting time is made shorter than 0-15 second. 

The feeding tube can only hold a few dozen balls, otherwise the pressure 
on the release mechanism gets too heavy. The electrically operated 
hopper shown can store many hundred balls and is very reliable. It 
consists of a hollow cone (its walls transparent for convenience) with 
a number of scoopers on its inside. The cone is slowly turned by a small 
motor and gear (4 revs./min.). Each scooper picks up some twenty balls, 
carries them to the top where it becomes suitably inclined because of the 
conical arrangement, and feeds them into the feeding tube. The latter 
has a crooked section which is hinged to the remainder by two short 
lengths of spiral spring (or soft rubber tube). When this section contains 
enough balls it presses upon a microswitch and turns off the motor 
which drives the hopper; when it is partly emptied the motor starts 
again. 

3. The propelling mechanism is a little billiard cue, driven electro- 
dynamically. The tip of the cue was made of polystyrene and seems 
quite satisfactory. The shaft is an aluminium rod, 2 mm. thick, threaded 
over its whole length (of about 7 cm.). It is clamped, with the help of 
four nuts, onto two V-shaped supports of beryllium copper foil so that 
it can move lengthwise only. Its other end is screwed on the driving 
coil (22 mm. diameter, 800 turns of 48 S.W.G. enamelled wire) which 
moves in the annular gap of a commercial loudspeaker magnet. 

When the cue strikes it gives the ball a velocity which is in a fixed 
proportion to its own (twice if the cue were much heavier than the ball, 
and perfectly elastic). The distance the ball will travel is proportional 
to the square of its speed (neglecting friction). In order to make that 
distance proportional to the pulse amplitude the speed of the ball (and 
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hence that of the cue as it strikes) must be proportional to the square 
root of the pulse amplitude. This will be so if the cue is uniformly 
accelerated from rest until it touches the ball, by a force proportional to 
the pulse amplitude. To achieve uniform acceleration, one has to avoid 
all damping and friction; yet the return motion ought to be damped 
.to avoid bouncing. This damping was produced by a 5 Kohm resistor 
parallel to the coil, the resistor being in series with a diode which becomes 
conducting during the return motion only. 
4. The valves of the electronic presorter are contained in one standard 
» chassis ; a second chassis contains a power supply of standard design, 
supplying +300 volts (up to 300 mA.) and —200 volts (up to 150 mA.), 
both stabilized. 


Fig. 3. 


Standard pulse generator and pulse lengthener. All resistances are in kilohm 
and capacities in p»F except where otherwise stated. M means megohms 
pF means puF; numbers given thus [20] indicate steady state voltages, 


4.1. (See fig. 3.) The input pulses must be positive and between 
7 and 50 volts in size. They come either from the input terminal or else 
from a multivibrator (valve V,) which produces standard pulses” of 
adjustable amplitude and recurrence frequency, for testing. After being 
fortified by a cathode follower (valve V,) those which are larger than 
a certain preset value (adjustable through R,) trigger the valve pair 3-4 
whereby valve 4 is rendered non-conducting for about 100 usec. At the 
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same time, as the potential of the control grid of V; rises with the pulse, 
the difference between it and the potential at A is amplified by V, and 
causes V, to produce current which tends to abolish this potential 
difference. Thus A follows the pulse until it reaches its maximum ; 
thereafter D, ceases to conduct, and C, remains charged to the amplitude 
of the input pulse until valve V, becomes again conducting. In this . 
way each input pulse (provided it is large enough) produces at the point A 
a voltage equal to its own peak voltage, remaining constant for about 
100 secs. and then returning to zero. Ry, is adjusted so as to exclude 


noise pulses, yet to include all pulses which we want to analyse. The °* 


Fig. 4. 
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Amplifier and biased amplifier. 


meter M, averages the current through V, and thus indicates the rate 
at which pulses pass this hurdle: this rate should not exceed about 
200 pulses per second or else some bias levels will shift. 

4.2. (See fig. 4.) Valves V, and V, form an amplifier with a gain of 
about 4 (adjustable through R,). The first half of valve V,, is a biassed 
cathode follower (which will not draw grid current, even for large pulses), 
feeding the second half, a triode amplifier of gain 5. The (negative) bias 
has coarse and fine adjustment (R, and R,) and is indicated on the 
meter M,. Pulses below a certain minimum size (which depends on the 
bias) are not amplified by this arrangement ; of larger pulses, the portion 
exceeding the minimum size is linearly amplified and appears at the 
point B, now as negative pulses. 


ee 
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4.3. (See fig. 5.) After passing through valve V,, which serves to 
adjust (through R;) the D.C. level, the pulses are supplied, in parallel, to 
the double triodes V,, and V,3. V,, forms a trigger circuit which will 
produce a certain negative potential on point D* if the potential on its 
left-hand grid drops below a certain level v,, adjustable through Rg. 
However, if the pulse is so large that it drops below v,, then V,, is 
triggered as well; the current in V,, is thereby turned off and point D* 
becomes positive again. In this case only a very brief negative pulse 
is produced, which is filtered out by the low-pass network R,C,R,C3. 
Thus a negative pulse of standard size, somewhat delayed and rounded 
off by the low-pass network, is produced at point D every time the 
input pulse has an amplitude between v, and v,. These “ control pulses ” 


Fig. 5. 


Discriminator. 


can be counted by an external scaler connected-to the terminal E. They 
serve, furthermore, to open a gate (see 4.5) which admits the pulse proper 
to the output valve. 

4.4. (See fig. 6.) The waiting time T is determined by the time constant 
(adjustable through R,) of the trigger circuit formed essentially by 
valve V,,, When triggered, this circuit renders conducting the right-hand 
half of valve V,;, thus:preventing further control pulses from getting 
through until the waiting time is over. If the control pulse does get 
through it is sharpened and reversed by the trigger cireuit Vigv lts 
leading edge places positive charge on the left-hand grid of valve V,,; 
its trailing edge triggers it and ‘“‘ shuts the door behind itself”, for a 
time sufficient to allow another ball to be placed in position to be fired. 
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Thus a negative control pulse arrives at point F only if the pulse at 
point B has an amplitude between v, and v2, and if the time elapsed 
since the arrival of the last control pulse is at leat T. Valve 18 acts as 
a gate which allows the pulse from point B to pass on only if a control 
pulse arrives at the same time at F. (C, is a neutralizing condenser.) 
4.5. (See fig. 7.) The output pulse, taken from the anodes of V,, and 
hence positive once more, now has to be lengthened to about 20 millisecs., 
longer than the longest time required by the billiard cue to complete its 
motion towards the ball. For this, the pulse is made to charge the 
condenser C; (by means of the cathode follower valve 19, and the 
diode D,a) while triggering the double triode V9, which discharges C; 


/ 
Time veto and gating circuits. 


on returning to its untriggered condition. The lengthened pulse is fed 
to the output valve V,,, while the (unlengthened) pulse triggers Vos, 
producing a standard square pulse which, through V,3, operates the 
ball release. | 


/ 
OPERATION. 


In order to get reproducible results it is important to keep both the 
balls and the board very clean. The balls are washed periodically in 
acetone, while the board is wiped with dry cotton wool to remove any 
dust. The distance of the tip of the cue from the resting steel ball must 
be sufficient (about 3 mm.) for the ball to be projected into the furthest 

' groove without driving the output valve V,, into grid current. 
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A number of electronic adjustments have to be made before the kick- 
sorter can be used. R, is first adjusted so that an amplification factor 
of 4 is obtained between the input pulse and the pulse appearing at the 
input to the biased amplifier (V,,). M, will then indicate how much 
voltage the bias removes from the original pulse height. Next R, 
(section 4.3) is‘adjusted so that the range of pulses admitted corresponds 
to 5 volts at the input terminals. When these adjustments have been 
made it is necessary to arrange for the smallest pulse to fire its balls into 
groove | and the largest one into groove 30. This requires a combination 
of adjustments of R; and Ri). R,; adjusts the ratio of the smallest to 


Fig. 7. 


Final pulse lengthener and output valves. 


largest pulse passed, R,) just varies the absolute magnitude of the output 
pulse. R, and R,, have to be adjusted alternately by successive 
‘approximations. 

When clean the kick-sorter will remain stable to one channel after 
the initial warming up period and will fire between 80 per cent and 
90 per cent of the balls into one groove with a standard pulse, the 
remaining balls being distributed over the two adjacent grooves. — sf 

One of us (S. G. F. F.) wishes to thank the Department of Scientific 
and Industrial Research for a maintenance grant. 
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LXIV. A Localizing Geiger Counter. 
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SUMMARY. 


The constant speed with which the discharge spreads along the wire 
in a Geiger counter is used to localize ionizing particles. Signals are 
obtained from the wire and also from two electrodes near its ends ; 
from these signals a pulse is derived whose amplitude indicates the 
position of the initial ionizing event. Provision is made for stabilizing 
the speed of spread. 


Tus instrument determines the position of the initial ionizing event 
along the length of a Geiger counter, with an error of less than 1 cm. 
Its operation is based on the fact that the discharge in a Geiger counter 
spreads along the central wire in both directions with a constant velocity 
of about 10 cm. per micro-second. The time which elapses between the 
start of the discharge and its arrival at one end of the counter is therefore 
proportional to the distance of the original ionizing event from that end 
of the counter. Furthermore the sum of the times taken by the discharge 
to reach both ends of the counter should be independent of the position 
of the ionizing event, and this serves as a check on the correct operation 
of the counter. 

To obtain all the required information about the discharge, two small 
auxiliary cathodes of about 1 cm. length and 2 mm. diameter are placed 
at the two ends of the anode. Their potential is adjusted to that of the 
equipotential which they occupy in the counter. At the start of the 
discharge a sharp negative pulse is obtained from the anode by 
differentiation, using a 1000 ohm resistance earthed through a 0-25 uF. 
condenser. Two positive pulses are obtained directly from the two 
auxiliary cathodes when the discharge reaches the two ends of the 
counter. All three pulses are fed into head amplifiers, which produce 
positive output pulses. These are then amplified by identical wide-band 
amplifiers which operate three fast trigger circuits (see fig. 1). The 
trigger circuit (6) actuated by the anode produces a negative square wave 
lasting about 20 microseconds. The trigger circuits (a) and (c) actuated 
by the auxiliary cathodes produce similar positive pulses lasting about 
twice as long. Pulses from the triggers (a) and (6) pass through an 
“interval converter’ which generates a signal p, whose amplitude is 
proportional to the time interval between the onset of the pulses from 
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(a) and (b) and hence to the distance of the ionizing event from the 
“upper” end of the counter. In the same way a pulse p, is obtained, 
proportional to the distance of the ionizing event from the other end of 
the counter. 

Fig. 2 shows the electronics in the blocks (a) and (1). The first pentode 


(head amplifier) is placed close to the counter ; the next three pentodes 


Fig. 1. 


Amplifier Interval 
and Trigger Converter 


Fig. 2. 


5t240 


Bs pF 


—180 
0 : 


form a fast amplifier which is coupled through a diode to the trigger 
circuit, comprising the next two pentodes. The last two pentodes and 
the last diode form the interval converter (1), which works as follows. 
The pulse from trigger (6) cuts off the valve (8); this allows the 
condenser C, to charge up through a diode and a 125 kw resistor until 
the positive pulse from (a) brings valve (7) into conduction. The voltage 
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in C, is kept small compared to 240 volts, and therefore rises linearly 
in the time between pulse (a) and pulse (b). Thereafter C, is discharged 
(much more slowly) through a 1 megohm resistor to the bias line. 

The electronics in block (c) is entirely similar ; in block (6) it differs 
in that the first pentode is used not as a cathode follower but as a phase 
inverter, with an anode resistor of 1-6 kw, and that the required positive 
signal for both interval converters comes from point A, not B, as shown 
wigs, 1. 

The propagation velocity of the discharge increases with increasing 
over-voltage. It is thus possible to control the over-voltage of the 
counter in such a way as to keep p,+p, constant. This stabilizes the 


Fig. 3. 


Pulse size distributions obtained from localizing counter. Diagrams 1, 2, 3: 


pulse Pr, with B-source at 5cem., 12cm. and 19cm. respectively, from 
lower ”’ end of counter. 


Diagrams 4, 5, 6: pulse $(p,+~.) for the same source positions as in 
diagrams 1, 2 and 3 respectively. 

Diagrams 7, 8,9: pulse p, with «-source at 2, 12 and 22 em. from end 
of counter. A bias of 6-9, 4-9 and 2-9 V. respectively was used to keep the 
group near the centre ; 1 volt bias means a shift of six channels. 


constant of proportionality relating both p, and p, to the position of 
the ionizing event, and also serves to keep the counter in the correct 
operating condition. To get the sum of the two pulses they are applied 
to the two ends of a resistor, from the centre of which a pulse equal to 
5(P,-+ Pg) is then obtained. The circuit used to control the over-voltage 
is similar to that designed by D. H. Wilkinson (1950) for the stabilization 
of proportional counters. 
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The counter used had brass walls with an internal diameter of about 
3-8 cm., and an active length of about 25 cm. ; it was fitted with a narrow 
mica window extending over the whole of its active length, had a tungsten 
anode of 0-1 mm. diameter, and was filled with 10 cm. pressure of argon 
and 1-5 cm. of alcohol vapour. To test it a beam of f-particles about 
2mm. wide was directed at the counter window at various points along 
the counter. The pulses were analyzed with a 30-channel kick-sorter 
(Frank, Frisch and Scarrott 1951) and the results are shown in fig. 3. 
The small subsidiary group shown on diagrams 2 and 3 is thought to be 
due to spurious counts originating at that point of the counter. 

The spread of the individual pulse groups is much greater than that 
of the groups due to p,+ , (see diagrams 4, 5 and 6); hence it could 
not have been due to any shortcomings of the instrument, but perhaps 
to scattering of the incoming f-particles at the window. To test this 
another similar counter was made with a thin enough window to admit 
«-particles which would be scattered very much less. Indeed the groups 
were very much sharper, as shown on diagrams 7, 8 and 9. In this case 
one channel width on the kick-sorter corresponds to about 8 mm. along 
the Geiger counter, and it is seen that localization is possible with a 
standard error of less than 1 cm. 

My thanks are due to Professor Frisch who suggested this work, to 
Mr. G. G. Scarrott for his advice about the electronics, and to the 
Department of Scientific and Industrial Research for a grant. 
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§1. INTRODUCTION. 


THE aim of the present treatment is to explain the development of 
deformation textures in polycrystalline aggregates from the behaviour of 
their single crystals, while at the same time attempting to overcome certain 
difficulties and limitations of previous methods. Its application to face- 
centred cubic metals led to the prediction of the polyerystalline stress- 
elongation curve, and, in a qualitative fashion, to the preferred orientations 
resulting from tension, compression, rolling, and wire drawing operations, 
(Calnan and Clews 1950, subsequently referred to as Part I.). It was 
shown that these predictions are in reasonable accord with experimental 
observations, although the latter are by no means comprehensive. 

It has been shown (Mises 1928, Taylor 1938) that, in the idealized case 
of homogeneous deformation, multiple slip, 7. e. slip on three or more 
systems, must take place in order to preserve the external form of the 
specimen and to maintain cohesion at the grain boundaries. It is now 
postulated that for this multiple slip to occur the stress system within a 
grain or part of a grain must be such as to give equal resolved shear 
stresses on all the operative slip systems. The simplest stress system satis- 
fying this condition is a single tensile stress, previously referred to in 
Part I. as the effective stress, T,, lying, in the case of a face-centred cubic 
metal, in a [100], [110] or [111] direction. In the initial condition with 
perfectly annealed material, at the instant of application of a small 
applied tensile stress, T,, the effective stress direction is coincident with 
the direction of this applied stress. As soon as elastic deformation takes 
place intergranular stresses generally come into play. If the applied stress 
is then increased and slip does not occur when the applied stress resolved on 
the most favourable slip system reaches the critical value for slip, then T, 
must have moved away from T, in such a direction as to reduce the resolved 
shear stress on this system. With further increase in the applied stress, 
provided no slip has occurred, T,, continues to move until it reaches one of 
the minimum positions of resolved shear stress, which are, of course, the 
points of multiple slip already referred to above. Having arrived at such 
a position, no further movement of T, is possible and there must be slip 
on all the symmetrically disposed equivalent systems when the critical 
stress is reached. Slip will frequently occur, however, before this con- 
dition is attained; there will be single slip if T, is within the reference 
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stereographic unit triangle, or duplex slip if T, is on a boundary between 
two equivalent slip systems. The operation of slip relieves the inter- 
granular stresses and the whole sequence above recommences. Thus the 
deformation of a single grain or part of a grain within the polycrystalline 
mass is visualized as a discontinuous process made up of elements of 
multiple, duplex, and single slip. It is clear that inhomogeneous deform- 
ation is implicit in this treatment. 

Due to the symmetrical disposition of the equivalent slip systems at the 
positions of multiple slip there is no grain rotation and the development of 
preferred orientation thus depends on the rotations resulting from single and 
duplex slip. Consequently, in order to produce the high degree of preferred 
orientation observed after extensive deformation, there must be a con- 
siderable amount of single and duplex slip, although some multiple slip is 
necessary to retain both cohesion at the boundaries and specimen shape. 

These concepts, described in greater detail in Part I., have now been 
applied to the more complex case of the deformation of a body-centred 
cubic metal where slip can occur on systems of the types {110}{111), 
{112}{111), or {123}<{111). Andrade and his co-workers (1940), who 
have made an extensive study of the slip systems in body-centred cubic 
metals, find that the slip direction 111) is invariate but that the slip 
plane depends on the temperature in such a way that it is the same for 
many metals for the same ratio of test temperature, T, to melting temper- 
ature, T,, (Table I.). 


me 


TABLE I, 


Slip Planes in Body-Centred Cubic Metals. 


Metal Slip plane AN 


Mo, W, Na 112 0-08 to 0-24 
B-brass, Mo, Na 110 0-26 to 0-50 
Na, K 123 0-80 to 0-87 


The systems operative in «-iron have been studied by Barrett, Ansel 
and Mehl (1937 a), who find that slip takes place on all three planes at 
temperatures between 20° C. and —195°C. Silicon ferrite, deformed at 
low temperatures or containing more than 4 per cent silicon, slips only on 
{110} planes. In conditions where all three planes are operative, 
however, the ratio of the critical shear stresses for {110}, {112}, and {123} 
slip, determined by Smoluchowski and Opinsky (1950), is 100 to 103 to 105. 

Deformation textures have been derived for each of the three slip planes 
operating independently, for the three systems simultaneously operative 
with equal critical shear stresses, and for the three systems with critical 
shear stresses in the ratios referred to above. 
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§ 2. Maximum RuEsoLveD SHEAR STRESS DIAGRAMS. 


Values of the function cos y cos A, where y and A are respectively the 
angles which the slip direction and the slip plane normal make with any 
particular direction of applied stress, have been tabulated at 1° intervals 
over ranges of «and f from 0° to 90° ; «is the angle of co-latitude and B the 
angle of longitude of the pole of the stress direction on a standard stereo- 
graphic projection. They have been calculated to three figures on 
Hollerith machines by the Mathematics Division of the National Physical 
Laboratory. 


Rigel 


Wen 


4Ro 480 WO 410 dbo $60 440 430 420 qotD 
110 
Maximum resolved shear stress contours, to three decimal places, plotted on 
a standard stereographic unit triangle for the systems {110}{111) and 
erates Some lines of quickest descent from the maximum at A are 
also shown. 


Contour diagrams of these functions for {110}, {112} and {123} systems 
have been published previously by Fahrenhorst and Schmid (1932). The 
justification for reproducing the diagrams based on these new tables lies 
not only in their greater accuracy but also in that they reveal some features 
overlooked in these earlier diagrams. 


{110}{111) and {111}{110) Systems. 


From the point of view of the evaluation of the function cos y cos A 
these two systems, are of course identical, although they are different 


ohio 
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crystallographically. The system ;111}{110), i. e. slip on {111} planes in 
(110) directions, has been considered in Part I. The related system 
{110} (111), slip on {110} planes in (111) directions, is required in the 
present discussion. The maximum resolved shear stress contours are 
shown in fig. 1. The most favourable system, that corresponding to these 
maximum resolved shear stress values, is the (101)[111] throughout the 
unit triangle. Along the edges of the unit triangle the maximum resolved 
shear stress is found on two symmetrically placed systems, while at the 
corners of the triangle the maximum value is on four, six or eight equivalent 
systems (‘Table IT.). 
Tasxe II. 


Most Favourable {110}<111) Slip Systems. 


Slip systems Number 


Unit triangle 

Boundary [100][110] 

Boundary [100][111] : ; 

Boundary [110][111] [111]. 

Point [110] [111], (011)[111], 


Point [111] 


Point [100] 


{112}<111) System. 

The maximum resolved shear stress contours for the {112}¢111) system 
are shown in fig. 2. This case is not as simple as the previous one, in that 
the same slip plane and direction are not the most favourable throughout 
the unit triangle. In the area [100]YZ the operative plane and direction 
are (211) and [111] respectively, while in the area [111][110]YZ they are 
(112) and [111]. Along the line YZ the resolved shear stresses for these 
two systems are equal. The most favourable slip planes and directions at 
triangle edges and corners are listed in Table III. It should be noted that 
the maximum number of most favourable equivalent systems is four, 
occurring at [100] and Y. 

The position of the boundary YZ between the two systems may be 
determined geometrically as follows. Consider the general case of two 
slip systems (p1917;)[P292"2] and (P3933) P4da74]- At any point of indices 
h, k, 1 the resolved shear stress on the system (p,%171)[ P2422] is equal to 
cos y cos A, where x and A are respectively the angles which the slip 
direction [p.qoo] and the slip plane normal [p,q,7,] make with the direction 
[hkl]. Cos x cos A is therefore equal to 


hp, +kathy hp stkqotlrs Bs 
PREGA (WERE) V/(B+ ath) (M+) 
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Fig. 2. 


50 ave BO CO 450 440 430 421 430 440-440 AG 41 
100 Y | a 


Maximum resolved shear stress contours for the system {112} 111), showing 
some lines of quickest descent from the maxima at B and C, and the 
boundary YZ between regions of different operative slip planes and 
directions. 


i 


TaB_eE III. . 
Most Favourable {112}{111) Slip Systems. 


Slip systems 


Area [L00]YZ 

Area [111][110]YZ 

Boundary [110]Z 

Boundary [110}[111] 

Boundary [1O0]Y 

Boundary [110]Y 

Boundary [111]Z 

Boundary YZ 

Point [110] 

Point Z 2TT)[111}, 
Point [111] ; » (211)[1TT)}. 
Point [100] y 117}, (211)[{117], 


1 
1 
1 
1 
2 
2 
2 
2 
2 
3 
3 
4 
4 


Point Y 111], (2U1)[117}, (112)[117), 


(egg eal 
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There is a similar expression, Q, for the resolved shear stress on the other 
system, and the boundary is given by P=Q, i. e. the equation 
(Rp P2—PsP4) +k*(Rqy42—9394) +2(Rryro—re 4) 
+hk(Rpide+Rpoqy —P314—P 44s) 
+h(Rpyro+Rper,—p3rs—P4"'s) 


+ kl(Rqyr2+ R¢or1—ds'4a—1ar"3)=0, foe UN ite lied oe (2) 
where R= (Gee ) 
(PI+Gi+11)(P3+43+73) 


This equation is homogeneous in the second degree and therefore is the 
equation of a cone. In the special case where the determinant 


A= (Rp \P2—PsP) | 2(Rpig2t+Rpq, (Rp ret+Rp,r, 
—Psd4—P 493) —Psl4—Pals) 
2(Rpiqet+Rpoy (Rq192—9394) 2(Rqyr2+Rqer, 
—Psd4—P Is) — sl 4—Ya"'s) 
2(Rp 72+ Rp, 2(Rayre+Rqor, (Rryr2—rsf4) 
—Psl4—Pa"'s) — 3" 4—a"'s) 


(3) 
vanishes the equation corresponds to a pair of planes. 
There is a further boundary corresponding to the equation 


ere Cee eee wen Werte; (4) 


which is also homogeneous in the second degree. The physical significance 
of this second boundary is that here the resolved shear stresses on the two 
systems are numerically equal but of opposite sign. The full expansion of 
(4) and its determinant differ from (2) and (3) only in that product terms 
with suffixes 3 and 4 are reversed insign. The intersections of the cones, or 
pairs of planes, with the edges [100][111], [100][110], and [110][111] of the 
reference unit triangle may be determined by substituting =k, /=0, and 
h=k respectively in (2) and (4). In the case of pairs of planes these 
intersections are sufficient to define the planes and hence the boundaries in 
the triangle. The form of the boundaries within the unit triangle in the 
ease of the cones may be determined by numerical solution in this region, 
and confirmed by the resolved shear stress tables. 
Applying these general results to the determination of the boundary YZ 
in the {112}{111) system, equation (2) takes the form 
h?—hk—2k?—31k+P=0, . . . .. . (5) 
which, since the determinant (3) does not vanish, represents a cone. The 
solutions of (5) for the conditions J=k, 1=0, and h=k are respectively, 
h—=2-56k or h=—1-56k, h=2k or h—=—k, and h=—1-781 or h=0-28l. 
A similar procedure for equation (4) shows that this boundary cone 
intersects I=k at h=0-59k and h—=—2-26k, 1=0 at h=0 and h=—k, and 
h=kat h=0-63l and h=—0-80l. These points of intersection are shown on 
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a standard stereographic projection in fig. 3. Thus the point Y is [210] 
and Z is [2-56, 1,1]. Fig. 2 has been drawn to be consistent with this 


information. 


{123}<111) System. 

The {123}111) system is more complex than either of those already 
discussed, since it involves three different most favourable slip systems 
within the unit triangle. The maximum resolved shear stress contours are 
shown in fig. 4 and the most favourable systems are listed in Table IV. 


Fig. 3. 


oOo! 


Standard stereographic projection showing the traces of the boundary cones 
corresponding to’ equation (2) for the systems (211){111] and (112)[117] 
equal in magnitude and sign (full line), and equation (4) for these systems 
equal in magnitude and opposite in sign (broken line) 


It is interesting to note that the corner points have the same number of 
equivalent systems as the {110}<111) case, and in addition there are six 
equivalent systems at X and four at Y. 

Consider now the boundary XY between the systems (312)[111] and 
(213)[111]. Equation (2) reduces to a pair of planes 


h+k—l=0, Pee. fe ee eee OG) 
h—2k—l=0, mah Sat ieee ieee eeermmat Ca) 
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while equation (4) also reduces to two planes 
h+k—l=0, ee een ee eee ee (7) 
h+l=0. re Ba eee) 


Tt will be seen that the two cones degenerate into three distinct planes 
Moreover, this degeneracy will always occur when the two slip peer 
under consideration have a common pole. The only plane which inter- 
sects the edges of the unit triangle is (6b), from which it may be deduced 


Fig. 4. 


Maximum resolved shear stress contours for the system {123}(111), showing 
some lines of quickest descent from maxima at D and E. 


that X is [311] and Y is [210], and further that the boundary is the great 
circle from [101] through these points to [111]. For the other pair of slip 
systems, (312)[111] and (312)[111], (2) and (4) become respectively 

3kl—hl=0, hl ee ee (5) 
and Sh? 2hk—R—2P=—0.  . .. « 6 » (9)- 
The coné (9) does not intersect the reference triangle edges and may be j 
neglected. (8) obviously represents a pair of planes 1=0 and h=3k, the 
traces of which are [100]Y and XW, where X is [811], Y is [210], and 
W is [310]. These boundaries are the parts of great circles. 


2U2 


624 BE. A. Calnan and C. J. B. Clews on the 


Three Active Systems. 


(a) Equal Critical Shear Stresses. The simplest case to consider when 
the systems {110}{111), {112}111), and {123}{111) are all possible slip 
systems is when the critical resolved shear stresses for the activation of 
these systems are equal. An indication of the most favourable systems 
may be obtained by the superposition of the shear stress contour diagrams 
for the individual systems, figs. 1, 2, and 4. The exact positions of the 
boundaries between the systems may be found by the method of the 
preceding sections. As may be expected, the system of boundaries within 
the unit triangle is relatively complicated (fig. 9). 


TABLE IV. 
Most Favourable {123}¢111) Slip Systems 


Slip systems Number 


Area [1l00]WX 
Area WXY 
| Area [110][111]XY 
| Boundary [100]Y 
| Boundary Y{110] é 
Boundary [110][111] 
| Boundary [111]X 
| Boundary X[100] 
| Boundary WX 
| Boundary XY y 
| Point Y 


| 


| — 
DOW NNWNWNN Ke = 


[111], 
[111], 


ns 


Point [110] 
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[111], (213)[111], (123)111), 
[111]. 

| Point X [111], 
[111]. 
Point [111] ub 
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Point [100] 
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2 
2 
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2 
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3 
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(b) Unequal Critical Shear Stresses. The preliminary results of 
Smoluchowski and Opinsky (1950) suggest that the critical shear stresses 
of the three systems are not equal. This case may be dealt with in the 
same straightforward manner with (2) and (4) modified to 


P= GO, ee ere 


where Cp and Cg are respectively the critical resolved shear stresses on the 
systems corresponding to the terms Pand Q. The changes in the positions 
of the boundaries for the observed critical shear stress ratios may be seen 
by comparison of figs. 9 and 11. 


ee 
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§ 3. TENSION AND COMPRESSION TEXTURES. 
Slip on {110}<111) System. 

Although the resolved shear stress contours are the same for the 
{110}<111) system and for the {111}{110) system discussed in Part L., 
the textures developed under tension and compression are interchanged. 
This is because under tension the stress direction tends to move towards 
the operative slip direction, the (111) in the {110}{111) case, and the 
{110} in the {111}<110) case, while under compression the stress direction 
moves towards the pole of the operative slip plane, the <110) and the 
{111} respectively. The directions of the rotations corresponding to 
single and duplex slip on the {110}{111) system in various parts of the 
unit triangle may be derived from a consideration of the operative slip 
systems listed in Table IT. These are shown diagrammatically in figs. 3 
and 6 of Part I. which now refer to compression and tension respectively. 
From these components, bearing in mind the relative probability of single 
and duplex slip, the general trends of rotation throughout the unit triangle 
are obtained. For example, at a point ¢ (fig. 1) it will be seen from the 
contours that the path of the effective stress, T.,, is first for a short distance 
to an adjacent point on the [100][111] boundary and thence for a much 
greater distance along the boundary to [111]. Thus it is reasonable to 
suppose that, in general, slip will occur more frequently when T, is 
situated on the boundary than when it is between the point ¢t and the 
boundary ; that is to say, duplex rotation is the more probable. For a 
point w in the area [100]AK near to A it will be deduced from similar 
considerations that single crystal rotation occurs more frequently. This 
leads to a tension texture as shown in fig. 8, Part I., and a compression 
texture as shown in fig. 5, Part I. Although minor components are 
predicted in these textures as a result of the sluggishness of rotations 
in the neighbourhood of the points of multiple slip, such components 
may well be absent after extremely severe deformation. 

The tension texture is principally a [110] texture with a spread between 
[311] and [110], which decreases with increasing deformation. There are 
no experimental data available for textures resulting form pure tension 
and discussion of this texture will be deferred. 

The compression texture is [111], [211] and [100], the last being the 
weakest component. Again it is probable that the [100] and [111] com- 
ponents decrease in importance with increasing deformation. 


Slip on the {112}<111) System. 

The single, duplex and multiple slip systems operative for the {112} (111) 
system are listed in Table ITI., and with the aid of these the general trend 
of crystal rotations in a specimen under tension may be derived (fig. 5). 
It is clear from this diagram that the texture will be mainly [110]. An 
interesting situation arises, however, in the region near to the minimum 
point, Z, in the resolved shear stress diagram. When the effective 
stress, T,, reaches Z slip occurs on three equivalent systems which are 
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not symmetrically disposed, resulting in rotation in the direction 
Z—+[1ll]. The part of this region lying in the triangle LCZ shows a 
complex behaviour, since the multiple and duplex rotations are in opposite 
directions quite distinct from the single crystal direction. Considering a 
single grain, or part of a grain, undergoing deformation, at some stages 
of the process it may suffer rotations corresponding to multiple slip and 
at other stages the reverse rotations corresponding to duplex slip. This 
clearly leads to a subsidiary texture which may persist until a very late 
stage in the deformation. Superimposed on these conflicting rotations 


Fig. 5. 


ur 


(ele) 
ite) 


(112}¢111 ) system. Tension rotations and texture. Regions of preferred 
orientation bounded by broken line. _ I, direction of single crystal rotation ; 
II, direction of duplex rotation ; m, direction of rotation due to multiple 
slip on unsymmetrically disposed systems. 


there is a general drift towards [111], the operative single crystal slip 
direction, which augments the [110] texture at the expense of this 
subsidiary texture. The point Y is also a minimum of resolved shear 
stress and a point at which multiple but unsymmetrically disposed slip 
occurs. Here, however, the multuple and duplex slip directions are the 
same, and the general trend of rotations towards the [110] as deduced 
from the duplex rotations is unaltered. The deformation texture is not 
modified therefore by the presence of this point of multiple slip. 
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The behaviour under compression is illustrated in fig. 6, and it may be 
seen that there is a double [111], [100] texture. In the main, orienta- 
tions within the area [111][110]YZ rotate towards [111] while those in 
the area [110]YZ move towards [100], so that the former texture 
predominates. The unsymmetrically disposed multiple slip systems at 
Y and Z modify this to a small extent, but do not give rise to subsidiary 
textures as in the case of tension. This is because, under compression, 
in the regions where the multiple and duplex rotations are opposed, 


Fig. 6. 


itt 


{112}{111) system. Compression rotations and texture. 


operation of a given rotation leads to a new orientation which favours 
further rotation of the same kind, whereas under tension operation of a 
given type of rotation leads to a new orientation which favours rotation 
in the reverse direction. 


Slip on the {123}<{111) System. 


From an examination of fig. 7 showing the grain rotations to be expected 
when slip takes place on the {123}{111) system under tension it is evident 
that the principal texture is [110]. The rotations near to EO EIANG 
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are complicated by the multiple slip at X, and they give rise to a subsidiary 
[311] texture which may be expected to persist until a late stage in the 
deformation process. The compression case is illustrated in fig. 8. Here 
there is a double [111], [100] texture, with [111] as the main component, 
and a subsidiary transitional texture around [311]. 


Fig. 7. 


All 


{123}<111) system. Tension rotations and texture. I, II, and m have the 
same significance as in previous figures. IL’ indicates a further possible 
duplex rotation. 


Slip on Three Systems. 


(a) Hqual Critical Shear Stresses. The rotations leading to the 
development of the principal tension texture, [110], and the subsidiary 
texture around [411] are shown in fig. 9. The subsidiary texture is due 
to the opposing influence of multiple and duplex rotations in this region. 
Since unsymmetrical multiple slip occurs at both the points U and V the 
subsidiary texture is rather more extensive than in previous cases. 
The compression texture is shown in fig. 10, with [111] as the main 
component, [100] as a minor component, and [411] as a transitional 
subsidiary component. 
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(b) Unequal Critical Shear Stresses. Comparison of figs. 9 and 11 
shows that the differences in rotations and the textures developed are 
relatively small despite the considerable changes in magnitude of the areas 
over which the individual slip systems operate. It appears that both 
the tension subsidiary texture and the compression [100] minor component 
are slightly reduced. 


he > 
Fig. 8. 


{123}<111}) system. Compression rotations and texture. 


§ 4, RoLLInG AND DRAWING TEXTURES. 


The rolling and drawing textures resulting from slip on any one or 
combination of the systems discussed may be derived in the manner 
described briefly in Part I. For example, with slip on three systems giving 
compression textures [111], [100], and [411], and tension textures [110] 
and [411], the [111] and [110] compression end points and the [110] 
tension end point may be simultaneously satisfied. The [41 1] transitional 
texture in compression satisfies both the tension end point [110] and, to 
a close approximation, the tension transitional texture at [411] as well 
as any spread from this point towards the [110]. It is therefore reasonable 
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to suppose that, although the [411] components are transitional both in 
tension and compression, they persist in the rolling texture which therefore 
includes, in order of decreasing magnitude, (411)[110 < 411], (111)[110], 


and (100)[110]. 


The rolling textures to be expected from slip on the various systems 


are listed in Table V. 


\0Q 


Three systems ; equal critical shear stresses. Tension rotations and textures. 
I, II, I’ and m have the same significance as in previous figures. m/’ is a 
further possible direction of rotation due to multiple slip on unsymmetrically 
disposed systems. 


System 


Fig. 9. 
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TABLE V. 
Rolling Textures. 


Textures in order of decreasing magnitude 


{110}<111) 
{112}¢111) 
{123} ¢111) 
Three systems 


[110 — 311}, (111)[110], (100)f110}. 
[110 — 211}, (100)[110]. 
)[110 <—311], (111)[110], (100)[ 110]. 
) {Lf 


l 
l 
l 
1)[110 <-411*], (111)[110], (100)f 110). 


(21 
(11 
(31 
(41 


* (411)[411] is not strictly a texture but (411)[511] is and the 
angular difference between the two is a matter of only about 3°. 
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The methods of Part I. may also be applied to wire drawing textures. 
These lead in all cases, after severe deformation, to a simple [110] fibre 
texture. 


§ 5. Discussion. 


The most striking feature of the tension and compression textures 
derived above is their similarity for all the slip systems considered ; it 
is only in the transitional components that they differ. 


Fig. 10. 
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Three systems ; equal critical shear stresses. Compression texture and 
operative slip systems. 


As for the face-centred cubic metals there are no data on pure tension 
textures. For compression, the experimental observations by Heger 
(1941) on an iron alloy containing 4-6 per cent silicon, which deforms on 
{110}111), showed a double fibre texture with [111] and [100] parallel 
to the axis of compression. The only other compression data available 
are for «-iron, where all three systems are active. Here again a double 
[111], [100] fibre texture was found with the [100] component the weaker 
(Barrett 1939). This double texture, with [100] the weaker, is Peed 
by the present treatment in all cases. It is expected, however, that 
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with only the {110}<111) system operative a [211] texture develops at 
the expense of the double texture after severe deformation. This has 
not yet been confirmed by experiment. 

In view of its industrial importance there is a considerable amount of 
experimental information on rolling textures in iron and steel, and a 
large measure of agreement from observer to observer. The most careful 
studies would appear to be those carried out in the Carnegie Institute of 


Fig. 11. 
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Three systems ; unequal critical shear stresses. Tension texture and 
operative slip systems. 


Technology (Gensamer and Mehl 1936; Barrett, Ansel and Mehl 1937 b; 
Barrett and Levenson 1941). They have examined the textures of rolled 
silicon iron alloy which deforms on {110}<111), and «-iron which deforms 
on the three systems. A typical (100) pole figure is reproduced in fig. 12. 
They concluded that there are two continuous sets of end points : (1) a set 
including the textures (100)[110], (511)[110], (311)[110], and (211)[110] 
indicated by the points A, B, C and D respectively, which have the 
[110] direction parallel to the rolling direction and the (100) plane at any 
angle up to 45° or 55° from the rolling plane ; and (2) a set of orientations 
including the textures (111)[110], E, and (111)[211], F, which have the 
(111) plane in the plane of the sheet and all possible positions of the rolling 
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direction in this plane. It was further found that, with {110}<111) 
slip or slip on all three systems, the only difference in the pole figures was 
that in the former case the scatter of (100) planes about the rolling axis 
was 55° to 60° compared with 40° to 45° for the three systems. It is 
obvious that no one of the predicted textures listed in Table V. can account 
for these continuous ranges of orientation. 


Fig. 12. 


RD 


(100) pole figure for rolled iron showing various ideal orientations : A, (100)[110]; 
B;(511)/ 110) 3¢C,(311)[110],¢ Dy (211)[110] s,.611.1)[110] ; F,.(111)[211). 
Heavy full lines represent a continuous set containing A to D. Circle in 
broken line is a second continuous set containing E and F. (Barrett and 
Levenson.) 


The introduction of a further simple assumption not only overcomes this 
difficulty but also explains the characteristically wavy slip lines in «-iron. 
It is that the ratio of critical shear stresses for the three slip systems 
varies from one part of a grain to another such that the polycrystalline 
aggregate may be considered as composed of units, some deforming on 
each one of the three systems and some on all three systems. From the 
considerable movement of the boundaries in figs. 9 and 11 it is evident 
that the variations of the ratios of critical shear stress from one region to 
another in a grain need not be large. The rolling texture of iron will thus 
include all those orientations listed in Table V., that is, those orientations 
illustrated by Barrett and Levenson (1941), with the trifling exception 
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that (411)[110] replaces (511)[110]. Moreover, the major textures of 
Table V. give the stronger set of components A to D, referred to in the 
preceding paragraph, with the exception of (100)[110] which is predicted 
as a weaker texture. This is probably a small point as the spread from 
the (411)[110]. on either side of the centre of the pole figure may well 
maintain this orientation. On the other hand it may be pointed out that 
Kurdjumow and Sachs (1930) did indeed find a saddle at the centre of 
their (100) pole figure. 

The careful experimental work of Gough (1928) showed that at low 
deformations the wavy slip bands in single crystals of iron can be resolved 
into short separate traces of slip planes of the now accepted forms. 
This he interpreted as duplex slip on the two planes nearest to the plane 
of highest resolved shear stress, but he was unable to explain why the 
deformation takes place by small elements of slip alternating from one 
plane to another rather than the production of two interlacing sets of 
straight slip bands. On the present assumption, however, the small 
segments of a wavy slip band correspond to parts of the crystal in which 
the different ratios of critical resolved shear stress produce different most 
favourable systems. ; 

The compression and tension textures predicted here now need re- 
consideration when applied to «-iron. The main textures are unaltered, 
but there are additional transitional subsidiary textures corresponding 
to slip on each of the three individual systems, 7. e. main tension texture 
[110] with subsidiary textures [211], [311] and [411], and main compression 
textures [111], [100] and subsidiary textures [211], [311], and [411]. 
These subsidiary textures are likely to coalesce into a general spread of 
orientation between [100] and [111]. 

The objection may be raised that in an iron alloy containing 4-6 per cent 
silicon the slip lines are usually very much straighter than in «-iron. 
Consequently this material must deform on a single system shown to be 
the {110}<111), and so could not attain the full rolling texture observed. 
However, Barrett, Ansel and Mehl (1937) have suggested that after 
severe deformation the same slip mechanism may operate in this alloy 
as in «-iron, which presumably implies that the effect of the silicon in solid 
solution in raising the critical shear stresses on the {112}{111) and the 
{123}<111) systems, although appreciable in the early stages, may be 
comparatively small after severe work-hardening with these as latent 
systems. The straight slip lines, of course, relate to observations on 
material after deformation slight in comparison with the rolling reductions 
under consideration. Furthermore, the predominance of the {110}<111) 
system may be expected to increase the component corresponding to this 
system in the final texture, the (211)[110 < 311] causing a greater spread 
of orientation along the equator of the (100) pole figure which has been 
confirmed experimentally. Thus it appears that the {110}{111) system 
plays a greater part in the deformation of this alloy especially in the 
earlier stages, and even in the later stages of deformation when the other 
systems are active in producing the full rolling texture. 
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Finally, it may be remarked that the predicted simple [110] fibre 
texture in cold-drawn wires is amply confirmed by experiments on all 
body-centred cubic metals so far examined. 


§ 6. CONCLUSIONS. 


The qualitative treatment of the development of deformation textures 
in face-centred cubic metals described in Part I. has been extended to 
those of body-centred cubic structure. The main features of tension, 
compression, rolling and drawing textures have been satisfactorily 
predicted. The detailed examination of the rolling textures of iron 
alloys presents some difficulty and necessitates the introduction ot the 
further assumption that in adjacent regions within a single metallic grain 
there are differences in the ratios of the critical shear stresses. This is 
justified in that it leads to the explanation of wavy slip lines in «-iron, a 
well established fact, and permits the prediction not only of the complete 
rolling textures of iron alloys but also some of the minor variations among 
them. 

It seems likely that the general concepts of this theory may be applicable 
to other aspects of the deformation of both single crystals and poly- 
crystalline aggregates. For example, there appears to be a close connec- 
tion between the alternative directions of rotation discussed here and the 
development of deformation bands. 

Further accurate experimental work on deformation textures is necessary 
before the. details of the theoretical predictions may be fully tested. 
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ABSTRACT. 


A description is given of the use of an organic scintillation counter in 
conjunction with a thin lens B-spectrometer to detect low intensity, high 
energy electron groups in the presence of a high background. Evidence 
for a 4-17 MeV. B-group in the radiations from Na” has been found, and a 
value of branching-ratio obtained. Measurements made on the secondary 
electron spectrum from Na+ show the presence of a high energy y-ray, 
possibly from the cross-over y-transition in Mg”. 


§ 1. INTRODUCTION. 


THe angular correlation measurements of Brady and Deutsch (1950) 
suggest that the spins to be associated with the 4:14 MeV., 1-38 MeV., and 
ground states of Mg” are 4, 2, 0, respectively, and that both the 2-76 and 
1-38 MeV. y-rays correspond to quadrupole transitions. This assignment, 
in the case of the 2:76 MeV. y-ray, is in disagreement with the-results of 
Rae (1949). who, from a measurement of the internal pair production 
coefficient, deduced that this was a dipole transition. Later measurements 
by Mims, Halban and Wilson (1950) on annihilation quanta arising from 
internally produced pairs gave results which were more compatible with a 
quadrupole transition. Another line of attack on this problem is the 
investigation of the cross-over y-transition by Bishop, Wilson and Halban 
(1950), who measured the energy spectrum of photo-protons in a deuterium 
filled ionization chamber counter. The crossover y-ray was not observed, 
and an upper limit of 310-7 has been placed on its occurrence. A 
turther investigation is that of Grant (1950) who attempted to measure 
higher energy f-rays corresponding to highly forbidden transitions to the 
ground and intermediate states using a thin lens spectrometer, with a 
coincidence counter pair for recording the focused f-particles. No high 
energy $-particles were observed and the upper limit placed on the inten- 
sity of the 4-14 MeV. f-transition was 1 in 104, 


§2. THE TECHNIQUE USED. 


The present work was pursued independently of that of Grant whose 
results only came to hand when the investigation was completed. 
Essentially, it was the same type of investigation as Grant’s except that an 
organic crystal scintillator was used to record the focused f-particles. We 
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have already given a very brief description of the technique (1950). The 
spectrum of pulses resulting from the irradiation of the crystal by the 
almost monokinetic electrons focused by the spectrometer is in the form 
of a symmetrical peak, with half width determined by the energy of the 
electrons and the crystal-photo-multiplier combination. In general this 
will be superimposed on a background due to cosmic radiation, and direct 
and scattered y-radiation. Scattered quanta and electrons will cause an 
increase in the pulse height distribution at small amplitudes. The pulse 
amplitude corresponding to the peak and the focused momentum will be 
related in the same way as energy and momentum of an electron. The use 
of the technique does enable one to pick out those pulses corresponding to 
electrons of the focused momentum from, in some instances, the much 
larger number of pulses arising from other causes, but recorded at the same 
value of coil current. 


§ 3. THe EXPERIMENTAL ARRANGEMENT. 


The thin lens spectrometer, adapted for ring focusing, gave a resolution 
of 5 per cent and transmission of 0-4 per cent when used with a 1 cm. 
diameter source. The sources were compressed blocks of sodium fluoride 
inside a thin aluminium covering, and with thicknesses ranging from 10 to 
40 mg./cm.”. The scintillator used was in the form of a clear block of 
napthalene +1 per cent anthracene 5 cm. in diameter and 3-8 cm. thick, 
grown by Mr. R. F. Jackson, Junior of A.E.R.E. and located outside the 
spectrometer behind a vacuum tight window of 7 mg./cm.? aluminium. 
The resolution obtained with this crystal was about 20 per cent for electrons 
of 1 MeV. The magnetic field due to the main coil was neutralized over 
the region of the photomultiplier by means of a set of coils carrying the 
main current, so that the gain was substantially independent of the 
focused energy. The response from the multiplier was amplified and 
analysed by means of a five channel Pulse Analyser Type 1074A. 


§4. THe Hieo EnERGy f$-SPECTRUM. 


Fig. 1 shows the pulse distribution obtained at values of momentum 
above that of photoelectrons from the 2-76 MeV. y-rays, and peaks due to 
the high energy f-group are clearly visible, superimposed on a back- 
ground due largely to direct and narrow angle scattered y-radiation. It 
was found that the pulse amplitude at the peak was directly proportional 
to the focused energy, as expected. Moreover a Fermi analysis of the 
results gave a straight line plot with an end point at 4-17 MeV. (fig. 2) in 
excellent agreement with the maximum energy of the f-transition to the 
intermediate state. Measurements extending over several half periods 
showed that the decay of the high energy f-transitions must be quite close 
to that of the main transition. An estimate of the relative intensities was 
made by comparing the slopes of the two Fermi distributions. This gave 
a value of 310-5 for the intensity of the 4-14 relative to the 1-38 MeV. 
B-rays, a result compatible with Grants. 
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The region beyond the end-point of the 4:14 MeV. f-ray was carefully. 
investigated but no trace of the 5-53 MeV. f-transition to the ground state 
was found. 7 values derived for the three transitions are 


fr(1-38 MeV.) =1-22 x 10°, 
fr(4:14 MeV.)=5-0 x 102, 
fr(5-53 MeV.) >1-2 x 1045, 


Fig. 1. 
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Variation of pulse analyser curves with focused momentum. 


The last two may be small by a factor of about two arising out of the 
assumption of a straight line Fermi plot for the forbidden transitions. The 
energy of the intermediate transition is in agreement with the order of 
emission of the two y-rays, as indicated by the results of inelastic scatter- 
ing experiments by Beghian, Grace, Preston and Halban (1950) and by 
Rhoderick (1949). 


§5. THe Hien Enerey y-Rapration. 
It was while performing a further check on the nature of the particles 
observed that it was noticed that the placing of a thick absorber over the 
source, so far from eliminating the peak in the pulse distribution due to 
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high energy £-particles, actually increased its height considerably. There 
can, we believe, be only one explanation to this phenomenon, and that is 
that the source emits high energy y-rays. A later check, however, using 
as a reference the Compton electrons due to the 2-76 MeV. y-ray ejected 
from the bare source, enabled us to show that in this latter case, at any 
value of the momentum, less than 10 per cent of the counting rate of high 
energy electrons could be ascribed to a high energy y-ray, and that the 
great majority were in fact due to the high energy f-transition. 
Measurements were continued with more intense sources, up to 10 mC., 
better resolution ~3-5 per cent and with saturation thicknesses of both 
lead and copper radiators. The shapes of the secondary electron spectra 


Fig. 2. 
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Fermi plot of high energy B-rays of Na. 


obtained (fig. 3) were quite different from the B-ray distribution, and had 
maxima in about the expected place. Some evidence of a photoline was 
obtained in the case of a lead radiator, but the value obtained for the energy 
of the y-ray from the photoline position, no correction being made for the 
finite thickness of the radiator, was some 8 per cent lower than the energy 
expected for the cross-over y-ray of 4:14 MeV. The value obtained from 
an extrapolation of the Compton distribution was £ per cent lower than 
4-14 MeV. The values obtained for the energy of the 2-76 MeV. y-ray 
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from the Compton distribution and photoline; however, were correct 
within 1 per cent after applying the 1 per cent correction in the energy of 
the photoline due to the saturation thickness of radiator, found to be 
appropriate at lower energies. It is felt, however, that the results are not 
necessarily incompatible with a y-ray of energy equal to the cross-over 
y-ray, since at 4 MeV., in lead, energy loss of fast electrons by bremsstrah- 
lung, is becoming important, and may effect the shape of both Compton 
and photoelectric distributions with a thick radiator. However, the 
possibility remains that the y-ray is some 4—7 per cent lower in energy than 
the cross-over, though this does not appear to be likely theoretically. 

By comparing the two Compton distributions the relative intensities of 
the high energy and the 2:76 MeV. y-rays was determined as 5x 10~, a 
result completely incompatible with that of the Oxford group. The results 
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of chemical analysis of both the sodium fluoride and a later source of 
anhydrous sodium carbonate, however, leave little doubt that the y-ray 
arises from the neutron irradiation of sodium, and measurements made with 
sources irradiated under cadmium, and in a lower neutron flux, exclude the 
possibility of a fast neutron reaction or neutron capture by Na” giving rise 
to the radiation. This, together with an approximate measurement of 
half-period makes it reasonably certain that the y-ray is associated with the 
decay of Na?4. 


§ 6. Discussion. 


Ifthe high energy y-ray and the 2-76 MeV. y-ray do in fact arise from the 
same level, the use of the Segre-Helmholtz formula, and the relative 
intensities, indicates that they differ by one order of multipolarity, and 
that consequently, the parities of the ground and intermediate states of 
Mg™ differ. This demands that the 1-38 MeV. y-ray be of odd multipole 


Highly Forbidden Transitions in the Decay of Na®4 641 


order, almost certainly dipole, and is incompatible with the angular 
correlation measurements. In the case, however, where the 2-76 MeV. 
y-ray is dipole, since it is known that the probabilities of emission of dipole 
and quadrupole radiation are about equal, it is possible that the cross-over 
y-ray could be octupole. This would require that the B-transition to the 
intermediate state, with a spin change of 2 and no parity change, should be 
second forbidden, and does not appear to be at all likely. The interpre- 
tation of the fr values of the various f-transitions is difficult because of 
uncertainty in the magnitude of the nuclear matrix elements. However, 
if the main f-transition is first forbidden, it is probable that the 4:14 MeV. 


f-ray is at least third forbidden, and the 5:53 MeV. f-ray at least fourth 
forbidden. 
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ABSTRACT. 


The thermal conductivities of German silver, stainless steel and 
constantan have been determined between 2 and 90° K. and the electronic 
and lattice components of the conductivity have been calculated. 

The electronic conductivity is limited by scattering by “ impurities ”. 
The lattice conductivity is limited by scattering of the lattice waves by 
electrons at the lower temperatures and by “impurities” at higher 
temperatures. 

A table is given showing the heat flow along specimens of each alloy 
resulting from temperature differences commonly met with in practice 
in low-temperature work. 


§1. INTRODUCTION. 


OF the various sets of experiments which have been made on the thermal 
conductivity of alloys at low temperatures the most detailed have been 
concerned with superconducting alloys. Measurements have been made 
on other alloys by Karweil and Schafer (1939), Allen and Mendoza (1948), 
Wilkinson and Wilks (1949), Schmeissner and Meissner (1950) and by 
Hulm (1951). 

The present experiments, on German silver, constantan and a stainless 
steel, were made originally because of their practical value and cover a 
wider temperature range than do any of the previous experiments. The 
thermal conductivity was measured between 2 and 90° K. and the 
electrical conductivity at a few temperatures in this interval. 

The electronic contribution to the heat conductivity has been estimated 
and, by subtracting this from the total conductivity, the lattice heat 


conductivity has been deduced. This is an appreciable fraction of the — 


total conductivity over the temperature range of these measurements. 
It is shown that up to about 25° K. the lattice conductivity is limited 
mainly by scattering of the lattice waves by electrons, but that above this 
temperature “impurity ” scattering predominates. 


§2. THe EXPERIMENTS. 
The apparatus used was the same as that described elsewhere (Berman 


1951), in which the temperature drop along the specimens is measured by 
two gas thermometers. The thermometers were fixed with Wood’s metal 
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to copper bands soldered to the specimens. The measurements of 
electrical conductivity were made in the same apparatus. 

The German silver was in the form of a }-inch rod, obtained from 
Messrs. Rollets, with composition Cu 47 per cent, Pb 2 per cent, Zn 41 per 
cent, Ni 9 per cent, and was measured as received. The grain size was 
estimated from a photomicrograph, using a grain contrast etch. The 
mean diameter was in the region of 0-02 mm. 

For the main experiments the length used was 5 cm., but two measure- 
ments were also made with a length of 2-5 cm. This was to verify (1) that 
the length to be assumed for calculating the conductivity is equal to the 
distance between the centres of the bands soldered to the specimen, and 
(2) that the corrections for heat losses are as small as had been calculated 
(23 per cent at 95° but less than $ per cent below 60° K.). The agreement 
between the two sets of experiments is very close (O and © in fig. 1) and 
it was not considered necessary to repeat this check for the other alloys. 

As the ratio of the room temperature electrical resistance to the residual 
resistance was very small (1-1) compared with the value found by Allen 
and Mendoza for a rod of similar composition, the specimen was given a 
similar heat treatment and heated to a dull red heat for a few minutes and 
allowed to cool slowly. The resistance ratio was, however, unaltered. 

The constantan specimen (Cu 60 per cent, Ni 40 per cent) consisted of 
317 enamelled wires of 36 gauge, in parallel. These measurements were 
made because constantan leads of about this gauge are frequently used in 
various low-temperature apparatus in which it is necessary to estimate the 
heat flow down them. The wires were soldered together at their ends and 
two copper collars were soldered round them at a distance apart of 3-5 cm. 

The stainless steel specimen consisted of three rods, of 2 mm. diameter, 
in parallel. This was an austenitic 18/8 titanium stabilized material of 
composition Ni 7-9 per cent, Cr 18-9 per cent, Si 0-7 per cent, Ti 1 per cent, 
C about 0-1 per cent. The austenite grains were estimated to be about 
0-01 mm. with some ferrite grains of about 0-002 mm. precipitated in the 
grain boundaries of the austenite. 


§ 3. RESULTS. 


The thermal conductivity and electrical resistivity of these alloys are 
shown in figs. 1 and 2. For the three alloys the thermal conductivity at 
the lowest temperatures varies as a power of the absolute temperature 
slightly greater than 1. At higher temperatures the conductivity becomes 
less temperature dependent. 

These results can be compared with those of other authors. The German 
silver specimen had a similar composition to both the German silver 
measured by Allen and Mendoza and the “ Silver bronze ” measured by 
Karweil and Schafer. Allen and Mendoza found, for a specimen with 
much lower electrical resistivity, that the thermal conductivity, measured 
between 1:3 and 4° K., was proportional to the cube of the absolute 
temperature. Hulm has raised some doubts as to the validity of the values 
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found by these authors. The present results agree well, however, with 
those of Karweil and Schafer who measured between 5 and 24° K. The 
extrapolated value at 113° K. is in good agreement with the value deter- 
mined by Lees (1908) on a specimen of somewhat different composition 
which had a simliar electrical resistivity. 

The values obtained on stainless steel are slightly higher than those of 
Wilkinson and Wilks between 10 and 20° K. on a specimen of similar 
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composition. They are, as would be expected, considerably lower than 
those of Karwell and Schifer who measured a steel with only 1 per cent 
of “ foreign ” atoms. 

It is of interest to note that below about 20° K. the thermal conductivity 
of this stainless steel is very similar to that of “‘contracid’’, measured by 
Karweil and Schafer, and of “chroman B2Mo’’, measured by Schmeissner 
and Meissner. The apparent rapid rise of conductivity above 70° K. 
found by the latter authors for chroman B2Mo has not been found for 
stainless steel. 


——— 
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There seem to be no previous measurements on constantan in the 
temperature range of the present experiments, but the values are similar 
to those found for another Cu-Ni alloy by Hulm, below 20° K, 


§ 4. Discussion or Resutrs. 


The heat conducted by a metal or alloy is the sum of the heat conducted 


by the electrons and by the lattice. The heat conductivity, K, can there- 
fore be written 


K=K,+K.,, 
Fig. 2. 
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where K, is the electronic heat conductivity and K, the lattice conduc- 
tivity. These two components are determined by the processes which 
scatter electrons and phonons (lattice waves) respectively and give rise to 
heat resistance. 
Heat conducted by the electrons. 

The electronic thermal resistivity, W,(=1/K,), is, to a first approximation, 
the sum of the “ ideal ’’ resistance W,, due to the thermal vibrations of the 
atoms, and the resistance due to the aperiodicity of the lattice field 
produced by the presence of “ impurities ”. According to Makinson (1 938) 
the ideal resistance of a metal at low temperatures should be proportional 
to the square of the absolute temperature. Although this ideal thermal 
resistance is not simply related to the ideal electrical resistivity, the 
thermal resistance which is due to impurities should be directly related to 


646 B. Berman on the Thermal 


the residual electrical resistance, pp, by the Wiedemann—Franz—Lorenz law 
and is therefore equal to p)/L)T, where Ly) is the Lorentz constant 
(4(ak/e)?=2-45 x 10-8 Watt-ohm/deg.’). From measurements by 


Fig. 3. 
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Griineisen and Goens (1927) of the thermal conductivity of copper and iron 
it is possible to estimate the order of magnitude of W; for the alloys 
measured. Using the extrapolated values of py. from fig 2, W, is about 
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I per cent of p/L)T for German silver and constantan at 90° K. and about 
_ 2 per cent for stainless steel. 

As W; decreases with decreasing temperature while Po/Liol increases, 
the electronic resistance can be ascribed to impurity scattering over 
the range of the present experiments. The electronic conductivity is 
therefore LyT/p9, and the lattice conductivity, K g must be K—L,T/pp. 


Heat conducted by the lattice. 


The lattice conductivity for the three alloys is shown in fig. 3. Up to 
about 20° K. it is proportional to T?, but above this temperature it 
increases more slowly and reaches a maximum value. In the case of 
German silver, the lattice conductivity decreases with increasing temper- 
ature at the highest temperatures at which measurements were made. 

In the same way as the electronic thermal resistance can be considered as 
the sum of the resistances due to different scattering processes, so the 
thermal resistivity of the lattice, W g 18 due to several processes which 
scatter the momentum of the phonons. The resistances due to these 
causes are not strictly additive but there are temperature regions in which 
one process is dominant and may be identified by the temperature 
variation of the resistance. 

Wy is the resistance due to ““ Umklapp ” processes by which the phonons 
mutually interact. At low temperatures Wy o T’e~*?7, where v is a 
constant and 6 is the Debye characteristic temperature (Peierls 1929). 

W>, is the resistance due to small scale lattice defects and impurities. 
At low temperatures W,, o T but is less temperature dependent at high 
temperatures. 

W,, is the resistance due to scattering of phonons at the grain boundaries 
and is proportional to T~? (Casimir 1938). 

These three scattering processes are also present in dielectric crystals, 
and the temperature dependence of the resistances due to them have been 
verified in other experiments (Berman 1951). In metals and alloys there 
is an additional resistance due to scattering of phonons by electrons ; 
this resistance may be written as W, and should be proportional to T~* 
at low tempe atures according to Makinson. 

For German silver the lattice resistivity, 1/K,, has been analysed into 

these various components, as shown in fig. 4. 
W, is estimated roughly from measurements on dielectrics having a 
similar value for 6 (~300° K.). Wy is of the order of 10 per cent of the 
total lattice resistance at 95° K. and this ratio falls off very rapidly with 
decreasing temperature. 

W, is calculated from Casimir’s formula 


Wz 4:5 x 10-662/rT? Watts-1 cm.-deg., 


where r is the linear dimension of the grains. Taking 0 as 300° K. and 
r as 0-02 mm., Wz is about 25 Watt-1! cm.-deg. at 2° K., which is only 
4 per cent of the total resistance at this temperature. At higher temper- 
atures this ratio falls off very rapidly. 
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In the temperature range of the experiments the lattice resistance can 
therefore be attributed mainly to Wp and Wy. At the lowest temperatures 
W,, is seen to be proportional to T™ ; this is the temperature dependence 
of the lattice resistance, Wy, due to scattering of phonons by electrons. 


Fig. 4. 
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§ 5. CONCLUSIONS. 


In the alloys, measured the heat conducted by the electrons and by the 
lattice are of the same order of magnitude. The electronic conductivity is 
limited by scattering due to the aperiodicity of the lattice, and the “ ideal ” 
resistance due to thermal vibrations is negligible below 100° K. The lattice 
conductivity appears to be limited mainly by scattering of the lattice 
waves by the electrons at the lower temperatures, and by the “ impurities ” 
at higher temperatures. 

In the case of German silver if the effect of the electrons were removed 
the heat conducted by the lattice would increase to a value which, below 
4° K., would be determined by the boundary scattering and not by the 
percentage of “‘ impurity ” in the alloy. 

Mendelssohn and Olsen (1950) have found that, in contrast to the 
behaviour of most metals, for some alloys of lead and for columbium the 
heat conductivity is greater in the superconducting state than in the 
normal state. To explain this they postulate a circulation heat flow 
analogous to that which has been used to account for the high heat 
conductivity of liquid helium II. Hulm (1950) has suggested that for 
these specimens the lattice conduction is appreciable-and increases when 
the specimen becomes superconducting owing to the reduction in the 
number of electrons capable of scattering the lattice waves. If this is so, 
then, although the electronic conductivity decreases, the total conductivity 
in the superconducting state may be higher than in the normal state. 

For the metals and alloys on which such experiments have been made, 
the impurity concentrations are much smaller than in German silver, so 
that the thermal resistance due to them will also be smaller. For relatively 
large grains (1/10 to 1 mm.) therefore, the lattice resistance due to even 
1 per cent of impurity would be less than the resistance due to boundary 
scattering at 2°K. If at this temperature the electrons no longer take 
part in both heat conduction and in scattering of lattice waves, then the 
total conductivity would increase by a large factor. 

It would therefore be of interest to determine the electronic and lattice 
conductivities, in the normal state, and the grain sizes of specimens which 
show the behaviour found by Mendelssohn and Olsen. It might then be 
possible to decide whether it is necessary to postulate a new heat transport 
mechanism to explain their results. 


ACP Te ND x. 


In the design of low temperature apparatus it is often necessary to 
calculate the heat flow, along tubes or wires, between parts which are at 
very different temperatures. If the temperatures at the end of a specimen 
1 cm. long are T, and T,, the heat flow per unit area per second is equal to 


Ts Tr I: ce 29 
K dT and may be considered as (T,—T,)K, where K is the “ mean 
at . . 

conductivity in the temperature interval T, to T,. 
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For convenience, the following table gives the heat flowing along 
specimens of German silver and stainless steel, 10 cm. long and 1 sq. mm. 
cross sectional area, and along a constantan wire of 36 gauge 
(2:9 10-2 mm.2), also 10 cm. long, when the temperatures at the ends 
vary from 80° to 0° K. 


Rate of Heat Flow along Specimen—in watts. 


Temperatures German Silver Stainless Steel Constantan 
at ends of 10 cm. long 1 mm.? | 10 cm. long | mm.? 36 gauge 
Specimen °K. cross-section cross-section. (2-9 x 10-7 mm.?) 
10 em. long 


8:0 10 30x10 ait x10 
4-8 x 10-3 1-7.4104 1-69<10-% 
aw al Ne 1 Lost ai xi) 
93x 105" Lite 3-5 x 10-8 
1-0x10-* 4:2 10-§ 3-5 x 10-7 
4:0 10-7 a sioe’ 10x10 


As these heat flows are often of interest in the determination of the 
volumes of liquid hydrogen and liquid helium which they cause to be 
evaporated, it should be noted that 1 watt evaporates about 100 c.c. of 
liquid hydrogen per hour and about 1 litre of liquid helium per hour. 
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LXVII. Observation of Cosmic Ray Events in Nuclear Emulsions 
Exposed in a Glacier at 3550 m. 


By J. B. Harpine*, 
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[Received March 7, 1951.] 


ABSTRACT. 


Photographic emulsions have been exposed at various depths under 
ice at 3550m. The variation in the numbers of slow p-mesons with 
depth is in fair agreement with a predicted variation. 

From the variation of the numbers of “ stars’ with depth, the best 
estimate for the absorption mean free path of the star producing 
radiation is 170--1) gm./cm.?. 

The variation of the numbers of 7+-mesons with depth is consistent 
with an absorption path length for nuclear interaction corresponding to 
a few times the geometric value. A large upward stream of 7-mesons 
has been observed and possible interpretations are discussed. 


§1. INTRODUCTION. 


In previous notes (Harding e¢ al. 1949 a and b) we have reported results 

on the absorption of the star producing radiation and on the variation 

of the number of z-mesons with depth in ice. In order to improve the 

statistics and the geometry, the experiment at the Jungfraujoch 3550 m. 

has been repeated. Experimental results for stars and mesons will be 
given. 

§ 2. EXPERIMENTAL DETAILS. 


In order to obtain good geometry a steel pipe (24-in. internal diameter, 
#-in. wall thickness) closed at the lower end was sunk vertically to a 
depth of 9m. into a glacier. Boxes of plates coated with Ilford G.5 
‘‘ Nuclear Research’ emulsion 200 thick were exposed vertically at 
depths down to 5 m. for 89 days. The pre- and post-exposure amounted 
to about four days spent mainly at sea level. Plates were also exposed 
at depths of 64 and 8}m., but these emulsions. were poured at the 
Jungfraujoch under a large thickness of rock. The pre- and post-exposure 
of these plates amounted to 2 hours at 3550m. The spaces between the 
boxes in the top 3 m. of pipe were filled with cylinders of paraffin wax, 
density 0-9 gm./c.c., which should produce nearly the same effects on 
cosmic ray particles as ice itself. 


* Present address, Atomic Energy Research Establishment, Harwell. 
+ Communicated by Sir George Thomson, F.R.S. 
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§ 3. EXPERIMENTAL RESULTS. 


In Table I. the frequency of occurrence of the various events are 
shown. We have assumed that the density of ice was uniform and equal 
to 0-9 gm./c.c. Only p-mesons of projected lengths greater than 130 u 
have been recorded; a p-meson being defined as one that does not 
produce a visible star nor that gives rise to a secondary meson. The 
numbers of 7*-mesons which stop in the emulsions and are travelling 
at angles 0-60°, 60-90°, 90-120°, 120-180° with the downward vertical 
at their points of entry, are given. All errors given represent standard 
deviation, based on the numbers of events recorded. Non-uniformity of 
the emulsion thickness may introduce additional errors of about 5 per 
cent. A random sample of plates were searched twice and by this means 
the overall efficiency for recording stars was estimated as 97 per cent, 
while for p-mesons the efficiency was about 80 per cent. 


Fig. 1. 


No. of p-mesons/10 degrees. 


Angle in degrees. 


Angular distribution of p-mesons. ‘The shaded areas must be added to the 
observed distribution (unshaded histogram) in order to correct for loss of 
tracks at large zenith angles caused by the thinness of the emulsion. 


The angular distribution of p-mesons arriving at the plate is shown in 
fig. 1. This distribution has been obtained by measuring 100 mesons at 
each of ‘the following depths :—0, 180 and 585 gm./cm.”. Since mesons 
making large zenith angles haveasmaller chance than those of small zenith 
angles of producing a track of at least 130, (and therefore of being 
recorded), this distribution is weighted to small angles. By considering 
the distribution in the angle of dip for various zenith angles approximate 
correction factors have been obtained and the corrected distribution is 
shown by the addition of shaded areas. The curve in this figure 
represents a cos? @ distribution per unit solid angle. 
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The p-mesons that we observe represent three types of events :— _ 


(a) 28 per cent of 7~-mesons leave no visible star on coming to rest 
in the emulsion (Bradner 1949) and these will be recorded as p-mesons. 

(b) Some will be »+-mesons produced from decay at rest of 7*-mesons 
in the immediate neighbourhood of the emulsion, the 7*-mesons having 
been locally created. | 

(c) The majority are the »+-mesons which form the hard component 
of cosmic rays. 

We shall now estimate the quantities (a) and (b) and obtain by 
subtraction those mesons which form the hard component. 

For every 100 7~-mesons which produce stars, 


28 (Gm ss) = 39 7~-mesons 
100—28 
stop in the emulsion and fail to produce stars. Only 80 per cent of 
the z--mesons which produce stars have ranges exceeding 130 and 
there is no reason to suppose that this ratio will be different for 
those 7~-mesons which produce no visible star. Because of searching 
inefficiency only 80'per cent of p-mesons of lengths greater than 130 u 
are recorded. Therefore for every 100 7~-mesons which produce stars 
39 x 0-8 x 0-8=25 7--mesons will be recorded as p-mesons. 

It is estimated that for every 100 7*-mesons stopping in the emulsion, 
37 .*+-mesons (of length greater than 130 ) will also stop in the emulsion 
due to the decay at rest of 7+-mesons in the surrounding material which 
is assumed to have the same stopping power as the emulsion itself. 
The slightly (~20 per cent) lower stopping power of glass, which, with 
the emulsion of the neighbouring plates, forms the bulk of the surrounding 
material, will only produce. a second order correction. From range 
measurements on the n+ of 7» events, it is estimated that we recorded 
only about 80 per cent of the 7+-mesons which stopped in the emulsion 
and as stated above, the searching efficiency for p-mesons of length 
greater than 130 u was alvo about 80 per cent. Thus for every 80 7*+-mesons 
observed to stop in the emulsion approximately 37 x 0-8=30 y.+-mesons 
are observed as p-mesons. 

Therefore, for every 180 7-mesons (80 7+- and 100 7~-mesons) observed 
to stop, 30-+25=55 mesons are recorded as p-mesons, and must be 
associated with the z-mesons rather than with the u+-meson of the hard 
component. By subtracting 55/180~30 per cent of the numbers of 
m-mesons at any depth from the numbers of p-mesons observed, we are 
left with the numbers of j.-mesons of the hard component. The values 
are shown in rows 16 and 17 of Table I. It is seen that only about 
4 per cent of the z-mesons travel upwards, and this upward stream can 
reasonably be attributed to Coulomb scattering. 

Working from the energy spectrum of 4-mesons at sea level, and the 
variation of the »-meson intensity with altitude, Sands (1950) has derived 
the expected form of the differential range distribution of j-mesons in 
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air at mountain altitudes. The experimental results for the numbers of 
#-mesons (Table I. row 16) travelling downwards and coming to rest in 
the emulsion, are plotted against depth of ice in fig. 2. Sands’ curve, in 
arbitrary units, has been corrected for the different medium (ice instead 
of air) and is shown in this figure. Thus assuming no p-mesons are 
produced in the ice, other than by the decay of z-mesons, the variation 
with depth in ice of the number of slow mesons, is in fair agreement with 
the predicted variation. 

As regards the absolute flux of .-mesons, Rossi (1948) has given the 
absolute value i, of the differential range spectrum at 10 gm./em.? at 
various atmospheric depths, 7. e. the number of .-mesons/day/steradian in 
the vertical direction, which stop in one gram of air at various depths in 


Fig. 2. 


Kep |. wd) 7 


ty =(-5 (after Rossi) 


n 
> 
2 
on 


No. of slow p-mesons/c.c/day. 


0 10C 200 300 400 300 600 700 800 


Depth in gm./cm.? ice. 


the atmosphere. At 3550 m. 7,=1-5 gm.~! diem.~1 sterad7?. This fixes 
the ordinate of Sands’ curve. In order to compare our intensity with 
Rossi’s value we have to correct for 

(a) the angular distribution of the mesons (cos? 6 per sterad), 
| (b) loss due to thin emulsions. We have recorded only mesons of 
length +130, whereas we require the numbers of mesons of length 
greater than zero, 

(c) relative stopping powers per gram of air and emulsion (1.gm. 
air = 1-7 gm. emulsion= 0-43 c.c. emulsion), 

(d) correction factor for searching efficiency is 1/0-8. 

We have recorded p-mesons arriving at the plate and coming from the 
upper hemisphere. It is to be expected that a photographic emulsion 
records the “ integrated intensity ” of slow mesons, 7. e. the number of 


2X2, 
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mesons recorded depends on the volume, but not the shape of the 
emulsion. Correction (a) is straightforward inasmuch as we have to 
compare the integrated intensity per sterad around the vertical with that 
of the upper hemisphere-correction factor is 


[27 cos? @sin @d6/["2rsin@d0 
limit aa 7S =Fz: 
e>0 ie 27 cos? 6 sin dé 2a 


Scattering, when once the mesons have entered the emulsion, will make 
their distribution more nearly isotropic and thus we obtain a factor to 
correct for (b) in the manner described by Lattes et al. (1947). The 
factor so obtained is 1-7 and we feel confident that this is correct to 
~20 per cent. We can now convert our intensities to Rossi’s units by 
multiplying by 3/27 x 1-70-43 x 1/0:8=0-44. The right-hand ordinate 
scale of fig. 2 is given in these units. It is seen that the experimental 
figures are approximately 25 per cent lower than those obtained from 
counter experiments, and the reason for such a discrepancy may be 
partly in correction factor (6) and partly in (a). Although the observed 
angular distribution of p-mesons more nearly obeys a cos?@ than a 
cos® @ distribution, part of these mesons are ~t-mesons formed in the 
immediate neighbourhood of the emulsion. These should be isotropic. 
Hence the p-mesons of the hard component sheuld be more strongly 
collimated. A cos? @ distribution, as has been suggested by Krauchaar 
(1949), would result in a correction factor for (a) of 4/27 instead of 3/27. 
Thus the discrepancy with counter experiments is not necessarily signifi- 
cant. From Sands’ curve, it would seem that the flux of fast u-mesons 
under 760 gm./cm.? of ice is ~1000/cem.?/day. 


Stars. 


Owing to the long exposure and to low development, the recording of 
tracks of particles at minimum ionization was not reliable. This 
introduces a small uncertainty, which should be the same at all depths, 
in the classification of a star. We have recorded stars which have (a) at 
least one track longer than 60, (6) at least three tracks with grain- 
densities greater than 30 grains/100 jy, 7. e. ~14 minimum, (c) a spur of 
less than 104 was not counted as a track. The number of stars per c.c. 
of emulsion per day are plotted on a logarithmic scale against depth 
in fig. 3. 

If the star producing radiation is coming vertically downward and is 
absorbed with a mean free path «, then the intensity at a depth a is 
given by 

Test Tere, os. Sieh ota dee ee eee en yy 

The results do not fit such a relationship too well, though the best 


value of « for the first 300 gm./em.? of ice would be 130-140 gm./em.2. 
There is a discrepancy with our earlier results (Harding et al. 1949 a) 
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from which a value of «200 gm./em.2 was obtained. The difference 
can only be explained by poor geometry, e. g. large holes in the ice in 
the first experiment. 

It may be argued that, because of an angular distribution of the 
star producing radiation, we should not expect exponential absorption. 
For an isotropic distribution at the ice surface of primaries, which are 
assumed to travel in straight lines, we would expect (see Rossi 1948) 


Pilg ee (ajax) Be( aay ee 2) 
or for a cos? @ distribution per unit solid angle 
T/Igp=e-*!*(1 —3(x/a) + $(x/x)?)—3(a/a.)®Ei(—aja). . . . (3) 
Fig. 3. 
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In order to obtain agreement between the experimental results and 
a theoretical curve over the first 300 gm./cm.? in the case of equation (2). 
There is then poor agreement for depths beyond 450 gm./cm.’. In the 
case of equation (3), if we take 160 <« <180 gm./cm.? good agreement is 
obtained between the theoretical curve and the experimental points 
even down to 765 gm./cm.? of ice. 

So far we have taken an exceptionally simple picture of the star 
production process inasmuch as we have assumed that all stars are 
produced by one component which is incident on the ice surface with 
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a certain (cos? 6) angular distribution. It seems likely that the star 
producing particles are not absorbed in a single nuclear encounter, or, 
what probably comes to the same thing, that in producing a star, the 
primary also produces one or more particles which can then act as primary 
particles for further stars, i.e. a cascade process. Provided these 
secondary particles have the same mean free path for star production 
and are produced at any depth with an angular distribution similar to 
that of the primary particles at that depth, then the above equations 
should equally apply for such a cascade process, but « no longer 
represents the absorption path length of the star producing radiations. 

It is thought that the nucleonic component is mainly responsible for 
the production of stars. However George ef al. (1950) have suggested 
a cross-section, o, per nucleon for star production by p-mesons of 
approximately 10-29 cm.2.. They have therefore shown ~0-1 stars/c.c./day 
should be produced in emulsion at 3550m. This of course is negligible 
compared with the 16 stars/c.c./day observed at the ice surface. Pro- 
duction of fast ~-mesons in the ice is negligible and hence at 760 gm./cm.? 
of ice there are ~1000 fast u-mesons crossing a square centimetre/day 
(see above). Thus the expected number of stars per c.c. per day produced 
by w-mesons is 1000Npo where N=Avagadro’s number and p=density 
of emulsion. This comes to 0-024 stars/c.c./day which is approximately 
one-third of the observed number at the lowest point. About 15 per cent 
of the stars at 585 gm./cm.? of ice could be formed by p-mesons. Thus 
even making allowance for stars formed by mesons, the experimental 
points still lie on a curve given by equation (3) with «=170+10 gm./cem.?, 

It is concluded that the experimental results are consistent with the 
assumption that the majority of stars are produced by primaries (nucleons) 
having an incident distribution cos? @ per steradian and an apparent 
absorption mean free path of 170-10 gm./em.?. 


a-Mesons. 


In a previous note (Harding and Perkins 1949 b) we obtained the 
variation in the number of 7-mesons with depth in ice, and concluded 
that they were produced with a mean range ~100 gm./em.? ice. Since 
then, an energy spectrum for the production of z-mesons in photographic 
emulsion has been obtained by Camerini et al. (1950). If we assume 
that this spectrum can be used to represent the production spectrum 
mm ice we can obtain certain limits for A, the absorption mean free path in 
nuclear interactions for fast 7-mesons in ice. 

It is assumed :— 


(a) the rate of production of mesons at a given depth is proportional 
to the star intensity at that depth, 


(6) mesons travel in straight lines, 


(c) the energy, or range distribution of mesons at production is 
independent of angle. Thus q(R) dR. f(@) dé is the probability that 
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a meson is produced with a range R to R+dR at 6 to 0-+-d@ with the 
downward vertical. The spectrum q(R) dR is obtained from Camerini’s 
spectrum using a theoretical range-energy relation. 


Let P(a, 0) dé 
aan o. of mesons stopping at depth (x, dx) travelling at angle (6, ay) 
No. of stars at depth (x, dz) 
where the star intensity (see above) can be approximated by 
I/Iy=e-2l, 


Let g(R) dR be the probability that a meson of range R gm./cm.2 ice 
will not decay in flight, 7. e. 


: fen Bike 
act le I, seam 


E being the kinetic, and E, the rest energy of a meson. The lifetime 7, 
is taken as 7 . 10~° secs. 


It can then be shown that 


[cos 6 Rcos 6 


P(x, 6) dd [rt (9)/cos 8 ee eB g(R)q(R) dR dé. 


Considering mesons in the angular range 0—60°, we can write as a sufficient 
approximation 


Rcosd 


ere. 
x)B « leks eB g(R)g(R) dR, 


where we take cos 6 as the average value of cos @ for mesons arriving 
at the plate (fig. 5). Cos 0=0-77. We take «=130 gm./cm.2. The form 
of P(x) can be obtained for various values of A, and these are plotted 
in fig. 4 together with the experimental results. We have fitted theoretical 
curves to obtain the best agreement with the experimental points by 
a method somewhat analogous to least squares. We choose the ordinate 
scale for the theoretical curves such that Ye?/o? is a minimum, where 


7 
e, is the difference between the experimental ith point and the curve, 
and o; is the uncertainty of the zth point. 

On this criterion, the experimental points best fit the curve for 
\=150 gm./em.2. The fit is about three times better than for A=75 and 
90 times better than for A=500gm./cm.2. We conclude that the 
variation with depth of the number of 7-mesons is consistent with an 
absorption path length for nuclear interaction corresponding to the 
geometric, or to a few times the geometric value of 70 gm. /cm.? (neglecting 
the effect of hydrogen). If the proper lifetime of 7-mesons is larger than 
we have assumed, say 10-8 instead of 7. 10~® seconds, then small values 
of X are still further favoured. 
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Assuming :— 

(a) mesons are produced only in oxygen in the ice ; 

(b) on the average the same number of mesons are produced in oxygen 
disintegrations, as in disintegrations of nuclei in photographic Se. : 

(c) the rate of nuclear disintegrations (and hence meson production) is 
proportional to the total geometrical cross-section of the nuclei, then 
we observe about 14 times as many mesons stopping as we would expect 
from a calculation in which a value of A=75 gm./cm.? was used. This 
shows that the numbers of mesons observed is of the magnitude expected : 
no particular significance is attached to the ratio of 1}. 

Since equilibrium between mesons and stars seems to be established for 
depths below 225 gm./cm.? ice, we should expect the angular distribution 
of 7-mesons to be constant at these depths. We are therefore justified 


Fig. 4. 


7-mesons stars x LO. 


Depth in gm./em.? ice. 


in classing together mesons at all these depths when considering the 
angular distribution of mesons arriving at the plate. This distribution 
is shown in fig. 5 in which a correction has been applied (shaded 
histogram) for loss of steeply dipping tracks. It is seen that approximately 
one-third of all 7-mesons were travelling in an upward direction when 
coming to rest in the emulsion. Unfortunately, complete data on the 
angular distribution at production of mesons have not been published. 
Brown et al. (1949) have shown that at 3550 m. in showers of more than 
five particles 93 per cent of the shower particles (the majority of which 
are mesons with energies above 60 MeV.) are directed downwards and 
nearly obey a cos @sin 6d6 distribution. We may assume the lower 


a 
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energy mesons to be isotropic. From an analysis of stars formed at 
70,000 ft. Perkins (private communication) concludes that in light 
elements about 10 per cent of shower particles travel backwards with 
respect to the incident star producing radiation. Due to the angular 
distribution about the zenith of the meson producing radiation (which 
is likely to be more collimated than the total star producing radiation 
since higher energy events are involved), the percentage of upward 
travelling shower particles at production will exceed 10 per cent. A liberal 
estimate of this percentage would be 20 per cent. We make the following 
assumption about the angular distribution of mesons at production and 
will then calculate the expected ratio of the numbers of mesons coming 
to rest after travelling in an upward cone of 60° to the number of those 

coming to rest after travelling in a downward cone of 60°. . 

Assume :— | 


(a) Mesons of energy below 60 MeV. are produced isotropically in the 
laboratory system. These form 12 per cent of the total number of mesons. 


Fig. 5. 
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Angular distribution of mesons. Shaded areas represent corrections for 
loss of steeply dipping tracks. 


(b) Twenty per cent of mesons of energy greater than 60 MeV. are 
produced isotropically in the upper hemisphere. 


_ (c) The remainder of the mesons are produced downward with a cos @ 
distribution per unit solid angle. 

Consider one meson produced in the upward 60° cone and one in the 
downward 60° cone due to group (a) mesons. Then it can be shown 
that there will be three mesons produced in the upward cone due to class (6) 
mesons and 18 in the downward cone due to class (c) mesons. With these 
assumptions we find that four mesons are produced in an upward cone 
to 19 in the downward cone. However, the downward moving mesons, 
being produced with higher average energies have longer ranges 
(electronic stopping) than the upward moving mesons. Therefore they 
should stand more chance of being lost through nuclear interactions 
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(for A=75 gm./em.2). Taking this into account, in a calculation somewhat 
similar to the above, it is deduced that one slow meson should be observed 
in the upward cone for every 3-6 in the downward cone. The experimental 
figure is 1 : 2-7--0-6. 

The calculations have been carried through on the assumption that 
mesons travel in straight lines. We have seen that the Coulomb 
scattering of -mesons is small (~4 per cent are back-scattered) and the 
scattering of 7-mesons is expected to be of the same order. We therefore 
conclude that for the given assumptions concerning the angular distri- 
bution at production of z-mesons, the observed ratio of the upward to 
downward flux of 7-mesons is not significantly different from the expected 
ratio. If, however, the results of further experiments should show that 
m-mesons are more strongly collimated than has been assumed above, 
then the observed upward flux would suggest that nuclear scattering of 
m-mesons is an important process. Indeed, preliminary results of 
Bernadini et al. 1950, suggest that 7-mesons are scattered in flight through 
angles of greater than 30°. with a cross-section of about one-quarter of 
the geometrical value. Scattering of this magnitude would lead to 
almost perfect agreement between the observed and ca_culated upward flux 
of z-mesons. 
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LXIX. CORRESPONDENCE. 


Meson Production in the Atmosphere. 


By J. D. Puttar and E. G. Dymonp, 
Department of Natural Philosophy, University of Edinburgh*. 


[Received March 12, 1951.] 


THE variation with altitude of the vertical intensity of the penetrating 
component of cosmic radiation has been measured in a series of balloon 
flights from Edinburgh, geomagnetic latitude 59° N. A telescope, ABC, 
of three counters with 8 cm. of lead absorber was used, giving a ground 
rate of 57-4+0-7 counts per hour. Radio transmission of the counts 
during flights and of the air pressure record, gave the necessary data. 
Particular attention has been paid to accuracy of pressure measurement, 
since the important information is obtained at the lowest pressures. 
A minimum pressure of 4mb. (maximum altitude of 124,000 ft.) was 
reached during these flights. 

Curve | shows the mean result of seven flights with identical equipments. 
The intensity, extrapolated to zero pressure, is 23 times the ground value. 
Pomerantz (1949) has published a similar curve, deduced from a single 
flight at A 52° N. and with 7:5cm. Pb. His curve, when normalized to 
ours at the ground, agrees within statistical error at the top of the atmo- 
sphere, but is some 25 per cent lower between 200 and 50 mb. At 27 mb. 
Schein, Jesse and Wollan (1941) found an increase of 16 times over ground 
rate, which is considerably lower than our result. Lack of published 
detail, however, prevents accurate comparison with our work. 

The counter arrangement ABC measures not only the flux of mesons 
and primary particles able to penetrate 8 cm. Pb., but also electron showers 
coming in from the side. The influence of these was measured in two 
further flights, by counting the fourfold coincidences between counters 
ABCD. The results are shown in Curve 2. At the lowest pressure of 
12 mb. the number of showers is rising very steeply and amounts to 20 per 
cent of the threefold rate ABC. 

Electronic showers arriving at large angles to the zenith are to be 

“expected at very low pressures. Such showers reach maximum develop- 
ment after traversing 100 gm./cm.2 of the atmosphere. Measured. 
horizontally the atmosphere has this thickness at a point where the 
pressure is 2-5 mb., thus the horizontal showers will increase in importance 
until this value of the pressure is reached. It is of course possible for 
the fourfold coincidence system ABCD to respond also to vertical showers 
containing at least one penetrating particle. However, the interaction 
length of protons is about 70 gm./em.?, so that a sharp rise in counting 
rate in the uppermost 20 gm./cm?. is not to be expected from this cause. 


* Communicated by Professor N. Feather, F.R.S. 
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We are therefore justified in regarding Curve 2 as due to side showers 
only, and by subtraction from Curve 1 we get the true vertical penetrating 
flux in Curve 3. This shows a maximum at 35 mb., representing ‘the 
production of mesons by particles of the primary beam. The agreement 
of our uncorrected curve with that of Pomerantz (who gave reasons for 
ignoring the shower correction) appears to be fortuitous. 

Extrapolation of Curve 3 back to zero pressure in order to determine the 
flux of primary penetrating particles remains doubtful until more accurate 
values of the lateral shower correction are available in the pressure range 
below 10 mb. Further work is in progress in this direction. 
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Curve 2: Coincidence rate of ABCD (lateral showers). 
Curve 3: Difference of 1 and 2, giving true vertical rate. 
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Cross-section for the Reaction Be (p, p2n)7Be at 156 MeV. 


By T. C. Ranpuz, J. M. Dickson and J. M. CASSELS, 
Atomic Energy Research Establishment, Harwell *. 


[Received March 22, 1951.] 


Tw the course of some experiments to measure neutron production cross- 
sections a ‘‘ sandwich” target containing two thick sheets of beryllium 
and three thin sheets of carbon was bombarded with 171 MeV. protons. 
The y-ray activities induced in both the target elements have since been 
followed for 104 days, using a Geiger counter and a standard geometrical 
arrangement. 

As usual (Dickson and Randle 1950, 1951), the radiation from the carbon 
part of the target quickly become characteristic of 7Be. The cross-section 
for the formation of this nucleus from carbon is 11-0-+1-0 x 10-27 em2. at 
the effective proton energy of 156 MeV. 

The activity of the beryllium part of the target could be analysed into 
two components, with half-lives of 10 days and 53 days respectively. 
The first component, which was much the weaker of the two, must have 
been due to an unidentified impurity. The second component was 
identified as 7Be, which has a half-life of 52-9 days. This assignment was 
confirmed by investigating the y-ray energy spectrum 57 days after 
bombardment, using a scintillation y-ray spectrometer (Hofstadter and 
McIntyre 1950). Mr. R. B. Owen, to whom we are indebted for this 
measurement, reported a single y-ray energy of 0-48 MeV. The y-rays 
emitted in the decay of 7Be have an energy of 0-485 MeV. 

A comparison of the activities of the two parts of the target showed that 
the ratio of the *Be (p, 2n) 7Be cross-section to the ™C (p, 3p3n) 7Be 
cross-section is 0-84-L0-04 at the effective proton energy of 156 MeV. 
The °Be (p, p2n) 7Be cross-section is therefore 9-2--0-9 x 10~®” em2, 

One of us (J.M.D.) intends to investigate the excitation curve for this 
reaction, which should prove useful for monitoring high energy proton 
bombardments of beryllium targets. 


We should like to thank Dr. A. Blainey for providing us with the 
beryllium targets. 
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The Reactions \B(ya)'Li and 1 B(yT')*Be. 


By Marcaret E. Catcrarr and E. W. TIrrEeRTon, 
Atomic Energy Research Establishment, Harwell f. 


[Received March 16, 1951.] 


As a result of experiments with the A.E.R.E. synchrotron, Goward, 
Titterton and Wilkins (1950) gave evidence for the reaction 


11B 4 ip 8H 9 He 11+ MeVir eee te 


and found that at least some of the events observed involved a two-stage 


process 
11B(/T)*Be* ; ®Be*->24He . . . . . . (2) 


the intermediate nucleus *Be being formed in the well known 3 MeV. state. 

Recent observation of the photodisintegration of Lithium-7 into a 
triton and an alpha-particle (Titterton 1950) suggests the alternative 
disintegration scheme 


11B(yo)"Li* ; "Lit>8H+4He. . . . ... .- (3) 


In the case where the “Li is left in its ground state a two-particle event 
would result according to the reaction 


1B hy 4He+‘Li—s-6 Mav. |... eer? 


instead of the three-particle stars of equations (1), (2) and (3). The 
reaction (4) has not been previously observed. 

To obtain information about these possible modes of disintegration 
Boron loaded E, emulsions have been exposed to Lithium gamma-rays 
from the 440 KeV. (py) resonance. Controlled under-development by 
the temperature method was employed to achieve discrimination between 
the distintegration products of reaction (4) if present. 

Although the experiment is still in progress and the statistics are poor 
certain conclusions emerge which are worth reporting at this stage. 

First, there is definite evidence for reaction (4) which is identified by 
the techniques discussed by Titterton (1950). A photomicrograph of 
one event is given in fig. 1 and the cross-section at 17-6 MeV. is estimated 
to be 2x10-*8cm.?. Twelve events have been observed, one of which 
is due to the 14:8 MeV. y-rays. 

Second, among the 750 three-particle stars which have been observed 
50 have been identified as resulting from reaction (1). The usual methods 
of the momentum and energy balance not only settle the assignment but 
enable a decision to be made as to which of the three particles is the 
triton. This is important in attempting to decide between reaction 
mechanisms indicated by (2) and (3). 


+ Communicated by the Authors. 
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Of the 50 events six result from the 14:8 MeV. and 44 from the 
17-6 MeV. y-ray. The cross-section for the reaction at 17-6 MeV. 
calculated in terms of the Cy3« cross-section given by Waffler and 
Younis (1949) is (0-6-+-0-3) x 10-28 em.2. 

Assuming all events to result according to the mechanism of equation (2) 
the Q of the break up of the recoiling *Be nucleus is calculated in each 
case and the result plotted in fig. 2. The histogram suggests two broad 


Fig. 1. 


Photomicrograph of event due to reaction 1B(y«)?Li. 
Observer : Miss M. Hart. 
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unresolved 8Be levels which can be identified as the well known 2-9 MeV. 
level and a level at about 4 MeV., evidence for which has been given by 
Green and Gibson (1949). There is no known level of *Be corresponding 
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to the small group at 5-5 MeV. In the case of the 2-9 and 4 MeV. levels 
the appropriate energy of a triton emitted at 90° to the incoming y-ray 
and the energy of the corresponding «-particle pair can be calculated 
from the conservation conditions for the 17-6 MeV. y-ray. The values 
obtained are indicated on the histograms given in figs. 3 and 4 which 


Fig. 3. 
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show the triton energy distribution and the energy distribution of the 
corresponding pairs of alpha-particles for events due to the 17-6 MeV. 
y-ray. In each case the observed distribution is consistent with the 
presence of events arising from the 2-9 and 4:05 MeV. levels. The groups 
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are not resolved but the distributions are consistent with the known 
resolution of the emulsion for low energy particles (Titterton and Brinkley 
1949) and the spread due to the angular distribution of the triton relative 
to the incoming y-ray, provided this is not markedly anisotropic. 

The 2-9 MeV. level of §Be is well known and has been much discussed : 
the 4 MeV. level was observed by Green and Gibson in experiments on 
the neutrons from the reaction’ Li(dn)*Be but they were unable to estimate 
the width of the state. However, as the only y-ray observed in the 
experiments of Bennett et al. (1941 and 1947) had an energy of 4:9 MeV. 
and the corresponding neutron group was observed by Green and Gibson 
it is reasonable to suppose that the 4 MeV. level leads to a pair of 
a-particles as indicated by the present experiment and has a width 
similar to the 2-9 MeV. level. The 4 MeV. level would then have even 
angular momentum or even parity or both. 

The mechanism indicated by equation (3) has been considered in a 
similar fashion. The known ‘Li level at 7-38 MeV. which breaks up into 
-a triton and an «-particle can be excluded immediately since none of the 
events show a low energy «-particle as would be necessary from the 
energetics of the reaction. The question then arises whether there is a 
lower level of 7Li which disintegrates by triton emission. The analysis 
is complicated by the fact that the first «-particle to be emitted is not 
known and each «-particle must be considered in turn. Two Q values 
relating to the break-up of *Li into a triton and an «-particle (second 
stage of reaction (3)) are found for each event one of which is correct and 
the other incorrect. A histogram similar to fig. 2 plotted with the results 
should show a peak superimposed on a background if a level of 7Li is 
involved. In fact, no such peak is observed and the experiment indicates 
that at least 85 per cent of the events proceed through the mechanism 
of equation (2). 

This result is consistent with recently published level schemes 
(Hornyak, Lauritsen, Morrison and Fowler 1950) for the mirror pair 
7Li and 7Be neither of which has a level between the low lying y-ray 
emitting states and the states at 7-38 and 7:02 MeV. respectively which 
are known to be unstable to «-emission. 


Our thanks are due to Miss E. M. Chellingworth and Miss M. Hart for 
microscope work and to Professor O. R. Frisch and Dr. D. H. Wilkinson 
for permission to use the Cavendish Laboratory H.T. set. 
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Observation of Growth of Cadmium Iodine from Aqueous Solution. 


By A. J. Forry, 
H. H. Wills Physics Laboratory, University of Bristol *. 


[Received April 23, 1951.] 
[Plate X XI.] 


AccorpING to the ideas of F. C. Frank (1949, 1950; see also N. F. Mott 
1949), the growth of crystals from vapour or solution can occur by two 
mechanisms. At high degrees of supersaturation (greater than about 
50 per cent), successive growth layers can be added to a perfect flat sur- 
face ; but for lower values, the free energy of surface nucleation is too 
high to permit this. At low supersaturation growth can only occur on 
surfaces where one of more screw dislocations terminate ; and the step- 
lines between growth layers should form spirals ending on the dislocations. 
Such patterns of monolayer step-lines have been observed by Griffin (1950) 
on prism faces of natural beryl crystals, using an optical microscope, and on 
paraffin crystals grown from solution by Dawson and Vand (1951), using 
the electron microscope. 

The appearance of certain faces of cadmium iodide crystals shown in the 
I.C.I. film on ‘* Crystal Growth ” (Bunn and Emmett 1949) suggests that, 
in addition to these monomolecular steps, steps of much greater height are 
formed. At present unpublished work by Verma and the author 
independently, suggests the same thing for carborundum ; these steps also 
form spirals. The spirals must end on a group of dislocations all of the 
same sign, the total strength or Burger’s vector of the groups being 
considerable. 

This letter records a series of photographs (Pl. X XI.) (taken with an 
optical metallurgical microscope) of the actual growth of a cadmium 
iodide crystal from aqueous solution, and shows the successive stages in 
the formation of the growth patterns of this type. 

Flakes of cadmium iodide were added to a drop of distilled water on a 
heated microscope slide until the solution was just under-saturated at 
60° C. Rapid crystallization around the edges of a cover slip placed over 
the drop formed a closed cell. The solution became slightly supersatur- 
ated as it cooled, and the growing crystal faces were examined in reflection 
with an ordinary metallurgical microscope. Step visibility was consider- 
ably improved by using a narrow illuminating pencil of light. 

A large number of small tabular crystals (about 100 microns across) were 
observed to grow quickly from the supersaturated solution, usually with 
the (0001) face parallel to the cover slip. Most of the faces exhibited 
patterns of diffuse concentric steps and the growth centres could not be 


* Communicated by the Author. 


A. J. FORTY | Phil. Mag. Ser. 7 Vol. 42, Pl. XXI. 
Fig. 1. Fig. 2. 


Figs. 1-5. (600 x) taken at 60 min. 
intervals, showing the upward 
growth of (0001) face of cadmium 
iodide crystal. 
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resolved clearly. In some regions of the solution, however, where the 
supersaturation was probably lower, giving a slower uniform growth, 
polyhedral layers were observed. The step-lines between these were 
sharp and sufficiently spaced for resolution of the growth centres, all of 
which showed patterns based on dislocation systems. A single screw 
dislocation dominating the growth of a crystal face produced the familar 
rotation spiral of step-line attached to it. A pair of screw dislocations of 
opposite hand, sufficiently close together to form a growing terrace, 
cooperated to produce a closed-loop pattern. Growth from two disloca- 
tions of the same hand gave a pattern of two parallel spirals, while growth 
from a system of three dislocations (two of one hand and one of the other) 
was indicated by a closed-loop pattern dominated by a single spiral. 
More complicated growth patterns have been observed from centres of 
many dislocations. On some very thin crystal plates, spirals have been 
observed on both surfaces with their origins one vertically above the other. 
This shows how a screw dislocation persists through the crystal layer 
lattice, producing a self-perpetuating step on both faces of the plate as it 
emerges. 

This “ spiral growth ” of the crystals is observed both for unpurified 
cadmium iodide and recrystallized cadmium iodide. In all cases, thin 
hexagonal or triangular plates grow rapidly from the supersaturated 
solution, and then a number of growth centres appear suddenly, to 
commence growth in thickness. Usually the growing face is dominated 
by one or two of the centres. 

It is likely that the lateral growth of the plates takes place by two- 
dimensional nucleation, from a high supersaturation, on the less closely 
packed (1010) of the hexagonal structure. This quickly lowers the 
effective supersaturation and further growth can only occur by the 
dislocation mechanism. F. C. Frank (1951) has suggested that the non- 
uniform distribution of impurities is sufficient to self-stress the thin plates 
beyond the theoretical yield stress. Buckling, followed by slip, produces 
the dislocation groups which have large screw components terminating in 
the (0001) face ; the crystal plate now grows in thickness from a low 
supersaturation. 

Photographs 1-5 are part of a sequence taken during the development 
of a particular crystal face. The step-lines have high visibility and show 
the trigonal character of the growth layers. The photographs illustrate 
the outward growth of the natural terrace between a pair of screw dis- 
locations dominating the face, until a closed-loop and a new terrace are 
generated. The cycle is repeated indefinitely, and the step-lines move’ 
outwards in closed-loops as the face develops upwards at the rate of one 
layer thickness per cycle. The growth of faces dominated by other 
dislocation systems has been recorded similarly. 

By illuminating the microscope with parallel monochromatic light, 
Fizeau fringes have been produced across this crystal face due to reflections 
from the glass-solution and solution-crystal surfaces. The thickness of 
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the growth layers has been calculated to be 750 A.+100 A. from the fringe 
displacement across the step-lines ; the steps are big and this accounts for 
their high visibility. 

As already indicated, this is not the growth of a crystal face from a pair 
of simple screw dislocations having unit Burger’s vector, generating 
unimolecular layers (these would have thickness 6-84 A. according to the 
X-ray data of Bozorth (1922)) as in the case of monolayer growth on beryl. 


The growth must be initiated from a pair of “ dislocation groups ”, each_ 


having a total Burger’s vector with a component normal to the face of about 
110 thicknesses of the minimum layer. The two groups cooperate to 
send out single growth layers 750 A.100 A. thick. Theoretically a 
dislocation group can exist of any total strength (7. e. any Burger’s vector 
normal to the face) and promote growth in layers of corresponding 
thickness. 

The two dislocation groups are not quite equal in strength and a small 
step corresponding to their difference extends outwards from the centre, 
crossing the main loops. At each crossing, the height of the main steps 
and consequently their rate of advance, changes abruptly, producing a 
characteristic sequence of pronounced kinks in the loops. Similar rows of 
kinks are observed in spirals on carborundum which grows in a like manner ; 
cadmium iodide can be used as a model for this type of growth. 

A further paper will give a more detailed account of these observations, 
with evidence that cadmium iodide crystals, like carborundum, can be 


polytypic. 
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